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VERTICAL EXTENT OF AURORAL ARCS AND BANDS! 


By B. W. Currie AND J. T. WEAVER 
ABSTRACT 


Measurements of the vertical extent of quiet arcs and bands from 181 auroral 
photographs, taken at Chesterfield (63.3° N., 90.7° W.), Coppermine (67.8° N., 
115.2° W.), and Saskatoon (52.1° N., 106.6° W.), indicate that these forms are 
confined to a narrow layer of the atmosphere immediately above the 100-km. 
level. The most frequent thickness of this layer is between 20 and 40 km. Values 
in excess of 50 km. rarely occur. The percentage of auroral time that the luminos- 
ity is restricted to this layer is uncertain. The ratio of the number of arcs and 
bands to the total number of observed forms suggests a value exceeding 50%. 
When the occurrence of each form is weighted according to the angular field 
covered by it, a value as high as 60% is indicated. 


INTRODUCTION 


Little is known about the range of heights within which each type of aurora 
commonly occurs, and the percentage of auroral time with each type as the 
dominant form. As a result, it is quite usual to assume vertical extents of 
several hundreds of kilometers in interpreting spectra, radio reflections, and 
magnetic and ionospheric disturbances observed during auroral displays. This 
is because heights for the upper limits of rays, as determined by parallactic 
methods, often range up to 300 and 400 km. and occasionally to over 1000 km. 
(5). 

A reliable determination of an auroral height from a pair of parallactic 
photographs depends on (a) well-defined edges to the structure and (6) an 
orientation of the edges roughly perpendicular to the direction of the base-line 
joining the photographic stations. The first condition is usually satisfied by the 
lower edges of arcs and bands, and the vertical edges of isolated rays. The 
brightness of arcs and bands rises rapidly to a maximum close to the lower 
edge and then drops slowly toward their upper edge until they are of the same 
general brightness as the sky. The brightness of rays is nearly uniform through- 
out their length. The second condition cannot be met by a single base-line for 
parallactic photography. Arcs and bands are nearly horizontal; rays are nearly 
vertical. A compromise is usually made, a direction for the base-line being 
selected that will give parallactic displacements making acute angles of about 
45° with the edges of both types of structures. As a result, a small error in 
locating the edge of aurora on one photograph relative to that on the other 
can lead to a large error in the parallactic displacement. 


1Manuscript received July 25, 1955. 
Contribution from the Physics Department, University of Saskatchewan, Saskatoon, Saskat- 
chewan. 
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An example will indicate the possible error in a measured height, due to an 
error Ap in a parallactic displacement p. For an auroral point 140 km. above 
the earth’s surface, and at an altitude of 20° in the most favorable orientation 
(the vertical plane through the base-line joining the photographic stations), 
p is about 1.5° for a base-line of 25 km. The relative error in a measured height 
is approximately Ap/sin p (see page 83 of Ref. 8). The scale of the projected 
drawings is usually one degree to a centimeter, so that a total error of as little 
as 0.5 cm. in locating the contours could result in a height error of roughly 
45 km. in this case. For the upper limit of arcs and bands the error in p is 
sometimes so large that computed heights less than the height for the lower 
limit of the same form are obtained. For the upper limit of rays the error is due 
largely to the small parallactic angles. 

A better approximation to the height of the upper edge of an arc is obtained 
by measuring the height of a point on the lower edge and the angular distance 
between this point and a point on the upper edge in the same azimuthal plane, 
and then by assuming that the two points are on an arc of a circle with radius 
equal to the distance to the lower point. If the vertical extent of the aurora 
of the afore-mentioned example were 40 km., the angular distance between 
bottom and top would be about 6°; and an error of 1 cm. in locating the upper 
edge on the projected drawing would result in an error of approximately 7 km. 
in its measured height. 

The assumption that elements of an auroral structure lie along arcs of circles 
in vertical planes is obviously faulty, and would lead to very large errors in the 
vertical extent of arcs for directions departing much from normal to them. 
Auroral rays lie along, or nearly along, the lines of the earth’s magnetic field, 
the exciting particles apparently spiralling down the magnetic lines. Presum- 
ably, the same excitation mechanism is applicable to homogeneous arcs and 
bands—suggesting that the vertical extent of these forms should be found by 
taking the difference between the upper and lower heights for points on a 
common line of the magnetic field. The fact that rays occasionally appear in 
such forms and are orientated in the direction of the magnetic field gives added 
support for this procedure. 


PROCEDURE 


Since a large number of photographs were available, both single- and double- 
station, that had been taken at Chesterfield (63.3° N., 90.7° W.), Coppermine 
(67.8° N., 115.2° W.), and Saskatoon (52.1° N., 106.6° W.), measurements 
were limited to quiet arcs and bands that crossed the geomagnetic meridian 
plane for each place. Evaluations for these cases were comparatively easy to 
make. They required measurement of the altitudes of the upper and lower 
auroral points in the plane of the geomagnetic meridian, a knowledge of the 
height of the lower auroral points, and a suitable chart on which the auroral 
points could be located and from which the vertical extent could be read. 
Much the same procedures were followed as for determining heights by double- 
station photographs (1, 5, 8). Outlines of the upper and lower contours of an 
arc or band were placed on tracing paper by projecting the photograph. The 
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tracing was placed on the appropriate ‘‘net”’ for azimuth—altitude measure- 
ments, and the altitudes of the upper and lower auroral points for the azimuth 
of the geomagnetic meridian were read from it. If the photograph was one 
from a set for parallactic measurements, the height of the lower edge was 
already known or could be determined; if not, a height of 100 km. was assumed. 
The reasons for selecting 100 km. for statistical studies are substantial, and 
have been discussed previously (2, 3). (Errors in the vertical extent due to 
height differences from the 100 km. value will be considered later.) A set of 
conveniently-spaced lines, representing the magnetic lines in the geomagnetic 
meridian plane for the place in question, was drawn on the chart generally used 
for determining the height and the distance to the horizontal projection of an 
auroral point (see Fig. 37 of Ref. 5). The magnetic lines were located according 
to the equation r/sin? a = constant, where r is the distance from the center of a 
magnetic dipole at the earth’s center to a point on the magnetic line and a is 
the angle between the straight line joining the afore-mentioned points and the 
axis of the dipole. Errors due to the fact that the field cannot be represented 
exactly by a dipole located at the center of the earth are small in comparison 
with errors from other sources. The magnetic line on which the structure was 
located could then be found from the height and the altitude of the lower 
auroral point. The height of the upper auroral point was found by following 
this line to the point with the altitude of the upper auroral point. The difference 
between the two heights gave the vertical extent. 

The percentage of auroral time that the luminosity is likely to be restricted 
to the range of heights given by the vertical extent for arcs and bands was 
found by counting the number of each form (arc, band, ray, drapery, etc.) as 
reported for Chesterfield during the Second International Polar-Year, 1932-33, 
and computing the percentage ratio of arcs and bands to the total number of 
forms. A similar count and computation was made for the 1932-33 Copper- 
mine data, although these were in less suitable form. Auroral glows were 
excluded, since they are of low intensity, and little is known about their heights. 

A more useful quantity is one indicating a combination of time and angular 
field covered by particular forms. The 1932-33 Chesterfield records were 
particularly suitable for estimating this quantity. Each observation included a 
sketch of the aurora on a zenithal-equidistant projection. A lucite disk of the 
same size as the projection was divided rather arbitrarily into 21 parts so that 
each part covered approximately the same angular field. The zenithal one was 
circular; the others were bounded by arcs of circles and radial lines. The disk 
was placed on each sketch, and a count made of the number of parts crossed 
or occupied by each form. This number was used to weight the occurrence. 
The percentage ratio of the weighted total of arcs and bands to the weighted 
total of all forms was then computed—again omitting auroral glows. 


RESULTS 


Figure 1 shows the statistical distribution of vertical extent for 144 quiet 
arcs and bands (52 Chesterfield, 33 Coppermine, and 59 Saskatoon) assuming 
a lower height of 100 km. Each point on the graph shows the number of occur- 











614 CANADIAN JOURNAL OF PHYSICS. VOL. 33 


45-49 


40-44 


35-39 


30-34 


25-29 


VERTICAL EXTENT 


20-24 
15-19 
10-14 


oe 5 10 15 20 25 30 


FREQUENCY 


Fic. 1. The frequency distribution of vertical extents of quiet auroral arcs and bands. 
Curve (a)fis from 144 single-station photographs assuming a lower height of 100 km.; curve (0) 
is from 37 sets of parallactic photographs with known lower heights. 


rences within the indicated range. On the same figure is shown the correspond- 
ing distribution for 37 arcs and bands at Chesterfield for which the lower 
heights had been measured by parallactic methods. 

While the number of these is small, the distribution is essentially the 
same as for the ones with the assumed height of 100 km. The lower heights 
for these ranged from 82 to 117 km., and the average value was about 102 km. 
The additional 2 km. in the average height explains in part the larger vertical 
extent for the maximum of this distribution. 

The percentage time ratio of arcs and bands to the total number of forms was 
54% at Coppermine, and 76% at Chesterfield. The percentage time ratio of 
arcs and bands to the total number of forms after weighting according to the 
angular field covered by each group was 67% at Chesterfield. 

DISCUSSION 

Several sources of error enter into the major distribution of Fig. 1, none of 

which are likely to alter the conclusion that the auroral luminosity for quiet 


forms is confined largely to a 50-km. layer of the atmosphere located mostly 
above the 100-km. level. These are (a) the difficulty in locating the upper edge 
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of the arcs and bands, (0) the possibility that the elements of an arc or band 
are not directed along the magnetic lines, (c) the assumption that the lower 
height was 100 km., and (d) the variation in density of the photographs. 

Numerous trials on different photographs indicate that error (a) seldom 
results in a departure of more than one degree in altitude from the average 
altitude for 10 trials on each photograph. Since altitudes both larger and 
smaller than the true altitudes are equally probable, all this error could do is to 
flatten out the distribution—leaving the maximum in the same place. 

The horizontal direction of arcs and quiet bands is nearly normal to the 
geomagnetic meridians (2, 3). Hence, measurements in the geomagnetic meri- 
dian plane would not have been altered appreciably by assuming other orien- 
tations for the elements, for example, along arcs of circles or vertically above 
the earth, as long as the angular altitude of the aurora was not large. For this 
reason, measurements were made mostly on arcs and bands within the altitude 
range from 15 to 25 degrees. 

While the effect of error (c) can be evaluated for particular cases, the total 
effect on the distribution of vertical extents is complex and can only be esti- 
mated qualitatively. Assuming that there is no error in the observed altitude 
of the lower edge, a difference between the true height and the assumed height 
of 100 km. affects primarily the distance to the selected point on the aurora. 
The difference between the true distance and the computed distance depends 
on the altitude of the aurora. The relative error in the distance, Ar/r, due to a 
departure AH km. from the assumed height of 100 km. is given by 


Ar AH 
r — -&(x—y) 
where x = (6272 sin’@+198.5)', 
y = 79.20 sin 0. 


6 is the angular altitude of the auroral point. The derivation of the above 
relationship is lengthy but not difficult. A constant radius of 6371 km. is 
assumed for the earth. Figure 2 shows the variation in the percentage error in r 
for altitudes up to 30° corresponding to AH = 5, 10, 15, and 20 km. For alti- 
tudes greater than 30°, the percentage error in 7 is practically the same as the 
percentage error in H. As the altitude decreases the percentage error in r 
decreases relative to that for H. A diagram showing several vertical extents 
for different values of r will show readily that the percentage error in the 
vertical extent is the same as the percentage error in 7 as long as the altitudes 
of the auroral points are the same. 

About one-third of many thousands of parallactic measurements on all 
types of aurora by Stormer (6) gave heights between 90 and 120 km. About 
two-thirds of the parallactic measurements at Chesterfield (1) gave heights 
between 90 and 130 km. for quiet arcs and bands. All measurements have 
shown a marked maximum of heights at about 105 km., and relatively few 
heights below 100 km. It is reasonable to assume that the true heights of the 
arcs and bands used in this investigation should have a similar distribution; 
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Fic. 2. Plots of the percentage error in the distance to an auroral point against its altitude 
for departures of 5, 10, 15, and 20 km. from an assumed height of 100 km. 


and that the most frequent vertical extent shown by Fig. 1 is not in error by 
more than 10%. It is likely that larger vertical extents than shown by Fig. 1 
occasionally occur, the assumed value of 100 km. biasing the distribution 
toward lower values. 

Since the auroral luminosity at the tops of arcs and most bands gradually 
merges into the general brightness of the sky, the density of the photographic 
image may affect the vertical extent, the denser photographs giving greater 
vertical extents. To test this point, the difference in density between the 
blackest part of each form and the general fog was measured for each photo- 
graph of the Saskatoon group. This was a homogeneous group with respect to 
photographic emulsion and developer. They were subdivided into four groups 
according to density and vertical extent, densities below the average for all 
forms and vertical extents below the maximum in their distribution curve 
being considered small. The number with small density and small vertical 
extent was 20; large density and small vertical extent, 9; small density and 
large vertical extent, 13; and large density and large vertical extent, 17. A 
chi-squared test of these values gives a probability slightly less than 0.05, 
suggesting that the larger vertical extents are associated with the denser 
photographs. An alternative hypothesis that the larger vertical extent occurs 
with the brighter aurora could not be checked. Neither exposure times nor 
visual estimates of the brightness of the aurora had been recorded in many 
cases. 
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Because of association of the errors (a), (6), and (c) with altitude of the 
aurora some significant variation of the vertical extent with altitude might be 
expected. A scatter diagram of vertical extent against altitude failed to show 
this. Chi-squared tests for various ranges of altitudes were also unsuccessful. 
Similar tests were made for variations of vertical extent with time of night, 
season, and latitude. All showed random variations relative to these variables, 
although a significant variation relative to the last was suggested by an 
average vertical extent of 27.9 km. for the Chesterfield and Coppermine 
photographs, and 31.0 km. for the Saskatoon photographs. 

Values greater than 50% for the percentage time ratio of arcs and bands to 
the total number of all forms are supported by similar estimates from data for 
other places. Fuller (4) in reports on auroral observations for the years 1931-34 
at College, Alaska, gives curves showing daily variations in the number of 
aurorae of each type seen at each half-hour. An average value of about 69% 
is obtained by measuring the areas below the curves for each type with a 
planimeter (glows excluded). Values given by Sverdrup (7) for the years 
1922-25 for the Arctic Ocean to the north of Siberia vary from 36 to 67% 
(glows again excluded). Conceivably, values of this quantity may vary 
according to the phase of sunspot cycle, and to the location relative to the zone 
of maximum auroral occurrence. While the idea of a percentage time ratio 
weighted according to the angular field covered by each form is not new, values 
for it do not appear in auroral literature. This may be due to insufficient 


observational data for most places. 
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ON THE EFFICIENCY OF AN ACOUSTIC LINE SOURCE 
WITH PROGRESSIVE PHASE SHIFT! 


By G. J. THIESSEN 


ABSTRACT 


The introduction of a small continuous phase shift along the length of a finite 
line source does not result in a decrease of energy radiated but simply causes the 
directionality pattern to change. When the phase shift is such that it corresponds 
to the phase shift with distance in a sound wave, then the directionality pattern 
has changed through 90° and the energy radiated begins to drop rapidly. At this 
point the slope of the energy vs. phase shift parameter curve increases with 
increasing frequency. 


I. INTRODUCTION 


Various special problems involving line sources have been discussed pre- 
viously in the literature. Olson (4) has considered directional characteristics 
of line microphones with a time delay introduced between adjacent lines and 
shows the effect on the directionality pattern. Others have considered the 
effect of a grounded slab (5) and a conducting half plane (6) on the radiation 
of a line source. Junger (2) has calculated the impedance ratio for various 
modes from an infinite cylindrical source in which the dynamic configuration 
of the cylindrical radiator is periodic in @ and z where 2 is measured along 
the axis. The subject of this discussion is the effect of a progressive phase shift 
along an acoustic line radiator of finite length on the directionality pattern 
and the efficiency. The term efficiency as used here is not the absolute efficiency 
but that relative to a line source which is in phase throughout its length. 


II. DIRECTIONALITY PATTERN 


Consider a line source of length / lying on the x-axis with the center at the 
origin as shown in Fig. 1. Let the source be such that the volume of fluid flow 
for unit length is given by 


Q(t) = Qoe**. 





Fic. 1. 


1Manuscript received July 21, 1956. 
Contribution from the Division of Applied Physics, National Research Council, Ottawa, 
Canada. Issued as N.R.C. No. 8782. 
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Then at large distances r at the point P the pressure due to an element dx 
of the source located at the origin will be given by (Ref. 3) 


os kpc ik(c t—r) 
dp = iz Qoe dx 


where k = w/c = 27/X, w is the angular frequency, \ is the wavelength, p is 
the density of the medium, c is the velocity of sound in the medium, and ¢ is 
the time. If the element is at x, then the path to P is effectively reduced to 
r—x sin 6. If in addition there is a phase shift along the source which varies 
with x so that the equivalent increase in path length is ax, then the pressure 
at P for an element at x will be 


. kpc ik(et—r+z8in0—ar) 
d =i{— e' ¢c r+z7sin dx. 
p 4nr Qo 


The total pressure p due to the whole line is obtained by a simple integration 
giving, after dropping the phase factor sin k(ct—r), 


__ p¢Qo sin[ki/2(sin @—a)] 
2ar sin 0—a ; 


[1 ] Pmax - 


The form of this expression is, of course, well known in optics and occurs as 
the directionality pattern of a long slit. The phase shift term a@ is less well 


a:0 A =0-5 


Oe 





Fic. 2. Directionality pattern for the sound pressure due to a line source as shown with a 
length 7 = \ and with a progressive phase retardation toward the right which is measured by 
the parameter a. 
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known and would occur if the slit were illuminated at an angle and even then 
the value of a would be limited. It is of interest to see the effect of varying a 
on the directionality pattern. In Fig. 2 are plotted such patterns for values of a 
indicated on each curve for the case where //A = 1.0. Although it is the 
pressure that is plotted here, it is obvious that the power radiated is not 
decreasing very rapidly until a = 1. 


III. TOTAL ENERGY RADIATED 


For the total energy E radiated we must integrate over all angles the 
expression 


Pinx/2pc, 


which gives 
pcQo™ +7?? cin? [k1/2(sin 0— -a)] 


mm 4n Jey (sin @—a) conte 
g TR(l—a) _- 2 
[2] e kipcQo sin'x oa 
8a —tR(l+a) X 


where x = k//2(sin@—a) and R = //\. Equation [2] can be changed into 
(see Reference 1, page 316ff.) 


becQ? Ft 5 


3 E= 
! 8 —rR(l+a) X 


where it is understood that J, is J, (x). For limits of zero and infinity this 
integral has a very simple form (7). However, the indefinite integral is some- 
what more complex. The general expression for integrals of this kind (1) 
contains an indeterminate term for this particular case, but can be evaluated 
from the derivative of the numerator and denominator leading to the final 
result 


aes ee 
16 3 FyJats, —*R(1+a) 


aes a ‘3 s (cd de i 


16 
oo j 4n—1 
~AR coe a) | 


Although equation [4] is a complete analytical expression for the total power 
radiated it is not an attractive one to use. However, its behavior as a function 
of a for the particular cases of //A = 5, 10, and 20 is shown in Fig. 3 where the 
energy radiated is compared with that for a = 0 and is plotted on the ordinate 
in decibels. The qualitative observation made previously that an appreciable 
decrease in energy radiated occurs only when a is near unity or greater is 
now more obvious. Furthermore, the slope of the curve at the point a = 1 
increases with increasing values of //X. 





[4] E = 





R(1—a) 


—7 R(1+a) 
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Fic. 3. The variation of the total energy output with a@ for //A = 5, 10, and 20. 


IV. CONCLUSIONS 


A continuous phase variation along the length of a line source affects mainly 
the directionality. Only when the rate of phase variation is comparable to, or 
greater than, that occurring along the sound wave (i.e. a phase variation of 2 
for a distance of \ along the line) will the energy radiated be appreciably 
decreased. Although the calculations were performed for an acoustical line 
source, very similar results apply to an electromagnetic radiator such as an 
antenna. 
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LOW ENERGY NEUTRON CROSS SECTION OF MANGANESE! 


By R. KrotKov? 
ABSTRACT 


The total neutron cross section of Mn® has three (predominantly scattering) 
resonances in the energy range 300-2500 ev. In the present paper these resonances 
are examined in detail from the point of view of the R matrix formalism. It is 
found that quite good (but not perfect) agreement with the experimental cross 
section curve is obtained by assigning (in order of increasing resonance energy) 
spins (2, 3, 3) and neutron widths (23, 11, 395 ev.) to the resonances. The addi- 
tion of a constant increment ARy = 0.16 X10-2 cm. representing the contribu- 
tion of distant, positive energy, J = 3 resonances to the R matrix element for 
s-neutron scattering improves agreement with experiment on the high energy 
side of the highest energy resonance. The (, y) cross section at 0.025 ev. is fitted 
well by a constant radiation width of 0.5 ev. However, with no set of parameters 
was it possible to fit both the experimental total cross section curve, and the 
observed coherent scattering cross section at 0 ev., although the sign of the 
latter was accounted for correctly. The discrepancy is 2.0 barns (calculated) com- 
pared to 1.7 barns (observed). 


INTRODUCTION 


From the point of view of neutron resonance reactions, manganese is a 
well-known and relatively simple nucleus. There is only one naturally occur- 
ring isotope, Mn®** (with 25 protons, 30 neutrons), with spin 5/2, and with a 
total neutron cross section (4, 19, 20, 14) as shown in the experimental curve 
of any of the figures.* The three peaks are at 337, 1080, and 2360 ev.; they are 
essentially scattering resonances (9, 1, 12), the only other process of importance 
being the (, y) reaction, (as is usual at these energies). Even this reaction is 
unlikely compared to elastic S wave scattering. The 337 ev. peak has been 
known for a long time (in 1947 Barbre and Goldhaber (1) estimated its width 
to be about 10 ev.), while the 1080 peak was found fairly recently (20, 4). The 
coherent and absorption cross sections at thermal energies are also known. 

In 1950, Harris, Hibdon, and Muehlhause (11) examined the then known 
experimental data to see if they could be fitted by a two-level Breit-Wigner 
formula (the 1080 ev. level was as yet undiscovered). They were able to obtain 
good agreement at several points, though they did not try to fit the curve in 
detail. They concluded that the 337 ev. peak had a neutron width of 21.6 ev., 
a radiation width of —0.6 ev., and a spin J = 3, while the corresponding values 
for the 2360 ev. peak were 304 ev., —0 ev., J = 2. 

A more detailed attempt to fit the experimental curve was carried out by 
Bollinger, Palmer, and Dahlberg (5), who took all three levels into account 
and found that best fit was obtained for a spin assignment J = 2, 3, 3, neutron 
widths I, = 22, 14, 336 ev., and a constant radiation width of about 0.6 ev. 
However, no choice of level widths was able to give complete quantitative 
agreement. 

The fast choppers at Brookhaven (20) and Argonne (4), together with work 
at Columbia (19) and Knolls (8), have recently given detailed data on the 


1Manuscript received May 31, 1955. 

Contribution from the Palmer Physical Laboratory, Princeton, New Jersey, U.S.A. 
2Present address: Palmer Physical Laboratory, Princeton, New Jersey, U.S.A. 
8This curve is taken from the compilation of Reference (17). 
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Mn* total cross section; the coherent and incoherent scattering cross sections 
at thermal energies are known from the early work of Fermi and Marshall (7), 
Bendt and Ruderman (2), and from the later work of Shull and Wollan (21). 
The activation cross section at 2200 meters/second has been measured by 
Pomerance et al. (15, 18), Colmer and Littler (6), Grimeland et al. (10), and 
Harris et al. (13). In the present paper an effort will be made to fit the 
experimental data from thermal energy to 4 kev. by the R matrix theory; 
the work of Harris et al. (11) and Bollinger et al. (5) has certainly shown 
that this theory can fit the curve in its essentials. It would be encouraging 
to verify that the full R matrix formalism can really “explain” the experi- 
mental data in detail. The next section will set up the equations to be used in 
numerical computations. 


MATHEMATICAL PRELIMINARIES 


As is well known, R matrix calculations apply to two-particle reactions and 
are carried out as follows (22). The total cross section is essentially specified 
by resonance energies £, and reduced widths 7, (the subscript s labels the 
channels, i.e. the possible two-particle states). These parameters determine the 
real, symmetric R matrix (with dimension equal to the number of channels) 


through the formula a 
ns Yt 
[1] Ry = ~ FE’ 
where E is the energy of the incident neutron (in the center of mass system) 
and R,, is the st matrix element of R. 
The next step is to find the scattering matrix U, which is given by 
1+iBRB—iC 
1—iBRB+iC° 
in terms of R and the diagonal matrices w, B, C. The latter are defined in terms 
of incoming waves J, in the various channels, all evaluated at the nuclear 


radius a,: 





[2] U=a 


WO, = 4 Nya 
[3] b, = [Im(1,/1,) 3, 
ay Re(J,/I,’) 
~ Im(/,/J. )- 


I,’ is the derivative of J, with respect to the channel radius 7;. w, B, and!C 


Cs 


also give the partial widths I, in terms of the reduced widths y;: 
ae 

4 Tye = ta ae 

[ ] As 1+c, Yas 


If R has more than two or three dimensions, U is difficult to calculate because 
it is hard to find the inverse of 1—7BRB+iC. However, once U is known, its 
elements give the experimentally measured cross sections directly: 

Tv - e e 
Os = 53 |Us:|" (s ¥ t, ie. a true reaction), 
8 


(5] 


Or = 7 |1—U,,|* (elastic scattering in channel s). 
& 


In [5], &, is the wave number of the incident neutron. 
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For us, there will be a channel for elastic neutron scattering (with / = 0), 
and a channel for each mode of y-decay (i.e. for each excited state of the 
residual nucleus), but these latter channels will be ‘‘weak’’ (matrix elements of 
R between these channels are small compared to the diagonal matrix element 
for the elastic scattering channel). Some modification of the standard scheme 
is required here, because, strictly speaking, the pair y-ray plus residual nucleus 
is not a two-particle system: it is not clear what should be used for the J and J’ 
(or B and C) corresponding to this channel. However, this difficulty can be 
avoided, and at the same time the inverse of the matrix 1—7BRB+i1C cal- 
culated, by proceeding as follows: 

The radiation reduced widths 7, (s now refers to a y-ray channel) are small 
compared to the neutron reduced widths y,1, which suggests that we break R 
up into two parts: 


R = Rote; 
0 | Rie Ris oo 
Ry | 0 Ro | Roe Res 
[6] Ro = » €= 7 Rg | Roe Ros 


| 
with e« small compared to Ro. We put s = 1 for the neutron channel. The 
matrices w, B, and C appearing in the expression [2] for U are unknown, 
except for the elements referring to channel 1. These elements are: 


o= e-e. 
[7] B = i, 
C=0 


(a being the radius of the Mn** compound nucleus). The elements 3,, c; 
referring to the y-ray channels we will have to determine somehow later; for the 
present we write: 


O ig 


U can then be calculated by using a power series expansion for the matrix 
(1—iBR)B+iC—iBeB)"! 
and keeping only the first order terms (in e). 
This expansion is (for two possibly non-commuting matrices A and in, 


n being small compared to A) 


(A —in)-! = A-\(1+iA-!+inA-inA=!+inA-inA-"nA-!+. ..). 
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In our case we put 
A = 1—iBR,»B+iC 
and n = BeB. 


For this A, A+A = 2 (A = complex conjugate of A), and the matrix of 
interest, w-!Uw!, becomes 


alia al ides i 
[9] aU ane = 4 +2iA "A T+... 
A straightforward calculation then gives the following expressions for the 
elastic scattering and (m, 7) cross section, os, and ony respectively: 
25 __ ,-2ika 1+7kRi | : 
[10a] Tue = 2 l—-e 1—ikRy 


[10d] = Liew = a 1+F? PE fa Ri 


s now refers to the y-ray channels only; the w, for the y-ray channels do not 
enter into [10] because they. are of absolute value 1 (cf. [5]). They would, 
however, enter into expressions for y—y scattering. 

From [10a] we see that o,, is just what we would get if there were no y-ray 
channels present, i.e. the weak y-ray channels introduce only higher order 
corrections to o-. The (m, y) cross section given by [100], on the other hand, 
is not yet in a form suitable for calculation. The right-hand side depends on the 
as yet unknown quantities },, c,, and on all the y,,. To see the dependence in 
detail, we can expand R,, 

Ri, = Ya1_Yas 


(11) Bate, 
interchange summations over s and X in [103] and so mee 


‘ _ <2 1 Ya Yul 

(12) om = "E TFPR EY a Eom Pe Sa nat 

From this equation we see that the vexatious b, and c, enter into o,, only in 
the form [b,2/(1+c,2)]? yas which from [4] is +(41T,,)?, the sign being that of 
Yas. Hence we can amalgamate [b,?/(1-++c,?)]# with y,, and consider the [,, as 
adjustable parameters with physical significance. However, op, still depends 
on all the [T,,; if the sum on s in [12] were independent of \ and yu (this sum 
could not then be a negative constant, because [12] came from [100], which is 
positive) we could treat it as a single-level width parameter and [12] would 
become quite a manageable expression. However, the sum is mot, in general, 
independent of \ and uw. We can set an upper limit on its absolute value by 
serene — 


| 3 Ss 
[13] | ) ce 2 cae YN Mas | < te ie. 2 Yas xX 14+¢, Yus 


- (8) (EY 





where T, = >, I, is the radiation width of level X. 





wa i is AR Ay LIM dn oe 





ine ls reins) 1 1s, 
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Since there are many y-ray channels into which the level \ can decay, we 
expect I’, to be independent of A, and this is indeed supported by experiment 
(16). Hence [13] becomes 


| b,” n\Vi/r,\!o rr 05 
it) Ter ore} B) a7 gre 


These considerations, together with the fact that capture processes are 
unimportant compared to scattering, suggest that the simplest way to make 
the right-hand side of [12] manageable is to replace the sum on s by a constant, 
3T',, and so get 


[15] i a aoe (x iM y 
' k 1+kRiv 2 y E,-E/° 

This is the expression for o,, that we shall use in the numerical calculations. 
I’, will be treated as a radiation level width whose exact relation to the 7, of 
[12] is vague but unimportant. 

In passing, we may note that [15] and [10a] are just the expressions (to first 
order approximation) one would get for on, and og, from a two-channel R 
matrix, if in the two-channel case we were to assume 72 (2 refers to the y-ray 
channel) to be independent of \, and put 


I, = b” 2 


2 ~T4+e2™ 

There is a difference between the two-channel and multi-channel cases, 
however. In the two-channel formula [15], o,, certainly goes to zero between 
resonances, while in the multi-channel case [100] it is very unlikely that o,, 
would ever become zero—for this could happen only if the R,, were all zero, 
and it would take a very special set of reduced neutron widths to accomplish 





this. 

The Mn data will, then, be fitted by formulas [10a] and [15]. Formula [10a] 
can be simplified by noting that in the energy range 0-4 kev., e?“** is practically 
1+2:7ka, and with this approximation [10a] becomes 
t(a—Ru)° 
1+k Riv * 


The actual observed total cross section is the sum of two cross sections, due 
to two R matrices—one, R- for J = 2, and another, R+ for J = 3. One of the 
set R+, R~ will be given by a two-level formula, the other by a single-level 
formula. As was first shown by Wick (23), the total scattering cross section 
will be an average of the scattering cross sections due to R+ and R-. At ther- 
mal energies the factor 1+7R),;? in [16] becomes 1 and we may write 
[17] Ose(thermal) = 42(g4(a—Rut)?+g_(a— Ri-)?) 


[16] Og = 40 


where g, = ;5, g- = are statistical weight factors. 

The coherent scattering between the two spin states (there is no isotopic 
incoherent scattering because Mn*5 has only one isotope) is described by the 
cross section 


[18] Ccon = 4n(g4(a—Riit)+g_(a—Ri17))?. 
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The sign of ¢eo, is defined to be the sign of a—g,Ri1t—g_Ri- (positive for 
hard sphere scattering). From [17] and [18] we find that the incoherent scatter- 


ing cross section is 
[19] Fine = Fse—Feon = 44 (g4g-)(Rut— Ri)’. 

We now proceed to the numerical calculations. With the seven parameters 
at our disposal (the three /,, the three reduced neutron widths 7, (ya: in the 


preceding formulas), and [',) we hope to fit the Mn curve from thermal 
energies to about 4 kev., and the following three data (taken from Ref. (17)): 


Osc (O ev.) = 2.0+0.1 barns, 
1.7+0.1(—) barns, 
13.4+0.3 barns. 


coh (0 ev.) 
Ony (2200 m./s.) 


SPIN ASSIGNMENT 


Of the seven parameters referred to in the last section, three can immedi- 
ately be found: the level positions E,. They are, of course, read off from the 
observed curve to be 337 ev., 1080 ev., and 2360 ev. 

A rough idea of the reduced neutron widths is found by fitting the curve in 
the neighborhood of each resonance. Thus at 300 ev., only the 337 ev. peak is 
important, so, taking a = 0.55X10-'? cm. (a reasonable value for Mn), 


2 2 
* ‘ eee {|| ae 
Srotar = 300 barns = 47 (0.55 nel 


or (y337)? = 290X10-!? ev-cm. 
Similarly, (v2360) = 2000 X 10—!2 ev-cm. 


The data in the neighborhood of the 1080 peak are poor, but a reasonable 
value for (y10s0)? is 100. With these approximate values we determine the spin 
assignment by reasoning as follows: 

(a) The 337 ev. and 2360 ev. resonances cannot have the same spin, for if 
they did, o,. (0 ev.) would be about 
290 2000 


a0. 2000) = 8.3 barns, 


o-(0) = 1 ae(0.55- 
where we approximate the spin factor (;5 or 4) by 3, (the effect of the 1080 ev. 
peak is negligible). This is much larger than the observed value of 2 barns. 

(6) The 1080 ev. and 337 ev. resonances cannot have the same spin, because 
if they did, their cross section would go to zero somewhere between them—at 
about 750 ev., say. At this point the total observed cross section would be due 
to the 2360 ev. resonance alone, and this is about 


i 
2 tx(0.55- anee = 3 barns, 


which is much less than the observed minimum of 16 barns. 

The spin assignment is hence either 2, 3, 3 or 3, 2, 2. 

To decide between (2, 3, 3) and (3, 2, 2) we appeal to the measurements of 
Bernstein e¢ a/. (3), who measured the capture cross section o,, of polarized 








EEE SOE NT SE ERA SOIRD I 
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neutrons on polarized Mn*5, and found that the antiparallel orientation 
(J = 2) had a little larger cross section than the parallel (J = 3). Now, at 
2200 m./s. ony due to the 337 ev. resonance is about 


1 4r 0.5 290 





2349x10* 2 (37)! ~ 11 bars, 
while o,, due to the other two is about 
1 4% 0.5 / 100 mae sas 
23.49X10° 2 (200. 2000 = 3.5 barns. 


We see that the 337 ev. peak predominates and so the spin assignment must be 
J = 2, 3, 3. This agrees with Bollinger et al. (5), but disagrees with Harris 
et al. (11). The spin assignment determines the peak heights of the various 
levels, the height of the J = 2 peak, for instance, being g_(44/k?) (k being the 
neutron wave number at resonance). The calculated peak heights are 3210, 
1410, 650 barns, the observed ones 2000, 240, 520 barns. This is reasonable 
agreement. (The second peak is so sharp that its full height is almost certainly 
not observed.) 
A FIRST ATTEMPT AT FITTING THE CURVE 

The spin assignment is now known (J = 2, 3, 3); the neutron reduced widths 
are about y,? = (290, 100, 2000) X 10-2 ev-cm. and I, is about 0.5 ev. (from 
a fit to on, = 13.4 at 0.025 ev.). Taking the nuclear radius a to is about 
0.55 X 10-!? cm., we can say the following about the curve: 

The (n, y) cross section can be neglected everywhere except at low energies 
(less than about 10 ev.). 

The 0 ev. scattering cross section is determined chiefly by the 337 ev. and 
2360 ev. peaks, and to a lesser extent by the 1080 ev. peak. o,,(0) is quite 
sensitive to the exact values of (7337)? and (y2360)? and, to a lesser extent, of a 
and (71080)?. The assignment y,? = 290, 100, 2000 gives too small a o,, at 0: 


_- .- 100 2000)" _ 
4r 12 (o.55- ay +4455 “(0 op ee = 1.6 barns, 


while the observed value is 2 barns. In addition, oj, is too small: 


35. (290_ 100 _ od _ 
Ome = Sag (20 1080 2360/7 ~ 9:02 barns, 


while it should be 0.3 barns. However, the sign of the coherent scattering is 
right: 
5 2¢ 2 
‘ine a 4 tf ee ome ‘ 


12 337 12\1080 * 2360 


As we move toward higher energies along the curve, but stay below the 
337 ev. peak, the 337 ev. and 2360 ev. peaks remain the most important ones 
in determining og, (and hence ojota; because o,, is small). As at 0 ev., the choice 
Yn? = (290, 100, 2360) X10-!? ev-cm. gives too low a total cross section. At 
100 ev. we calculate 3.8 barns against the observed 4.7 barns. 
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Moving again to higher energies, in the immediate neighborhood of the 
337 ev. peak our present choice of reduced widths, of course, gives rough 
agreement. Between the first and second peaks the J = 2 and J = 3 cross 
sections are about the same, and for the first time the 1080 ev. peak has become 
important. At 750 ev. we calculate ojota to be 15.6 barns, while the observed 
value is 16.0, so our choice of reduced widths gives about the right cross 
section here. 

The 1080 ev. peak is so sharp that the data in its neighborhood are not good; 
we shall not try to fit them any too well. Since the second and third peaks 
have the same J, 3, their cross section must be zero somewhere between them, 
and at this minimum the total cross section is the cross section due to the first 
peak only. For (7337)? = 290, o;ora: is about 4.2 barns at this minimum while the 
observed value is 7.5 barns. 

Passing on to the last peak, where, of course, the last level is dominant, 
(y2360)? = 2000 gives rough agreement with the experimental data, except on 
the high energy side (>3000 ev.). Here the calculated cross section is too high: 
at 4000 ev. we calculate 25 barns and observe 16.5. In this energy region the 
337 ev. peak contributes a roughly constant cross section of about 2 barns, 
while the 1080 ev. peak is quite negligible. 


ATTEMPTS TO IMPROVE AGREEMENT 


To obtain agreement at 0 ev., we must increase either (7337)?, (v2360)* 
(possibly (71080)?), or both, and also increase the difference Ri+—Riy~ at 0 
ev. The condition ¢,.(0) = 2.0, and ogo, (0) = 1.7(—) leads by straightforward 
algebra to the relation 

a—R,,*(0) — .242 /—.498 
or ( X< 10-** cm. 


a— Rj, (0) — .549 — 19] 


Taking a = 0.55 X10"? cm., this means (7337)? must be either 250 or 370. Both 
of these possibilities give a poor fit to the first peak, as Figs. 1 and 2 show. In 
these curves, (y1080)? and (y2360)” were determined by insisting that the mini- 
mum between the second and third peaks falls at the experimental minimum, 
1250 ev. Although when the curves were drawn this seemed justified, further 
reflection indicated that this was a mistake, because the experimental data 
around the minimum are not sufficiently reliable to warrant placing this much 
emphasis on them. I’, was determined by fitting ¢,, = 13.4 at 2200 m./s. 

It seems clear that in Fig. 1, (7337)? is too small, (y10s0)? too large, and 
(y2360)? about right, except for the 4000 ev. point. In Fig. 2, (7337)? is too 
large, (y10s0)? possibly a bit big, and (72360)? definitely small. 

In Fig. 1, the nuclear radius was taken to bea = 0.53, which does not make 
much difference to the curve above 200 ev., but improves the low energy fit. 
We could try to improve Fig. 1 by decreasing (71030)? to 170, which would give 
a fit at 750 ev., but also decrease o,,(0) to 1.5, which is too much of a drop, 
and at the same time leave the fit around 200 ev. poor. 
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Fics. 1 to 5. Comparison of the observed and calculated total neutron cross sections o; of 
Mn® as a function of neutron energy. Not all of the experimental points are shown; those 
appearing are intended only to give a general idea of how much they scatter. The triangles at 
the bottom give the experimental resolution. Using the parameters indicated, with a being the 
nuclear radius, the calculated cross sections at 0.025 ev. (neutron velocity 2200 m./s.) are as 
shown in the bottom left-hand corner. These are to be compared with the experimental values 
Ge = 2.0 barns, goon = 1.7 barns (—), ony (0.025 ev.) = 13.4 barns. 


In the next attempt we abandon our insistence that ogo, be 1.7(—), and pick 
(y337)? = 300, about right to fit the 337 ev. peak alone, but increase (72360)? 
to 2180 to get o,.(0) = 2.0 and then decrease (71080)? to get a fit at 750 ev. 
I, is determined as usual by fitting o,, = 13.4 barns at 2200 m./s. and ais taken 
to be 0.55X10-!2 cm. The results (Fig. 3) are not bad, but ogo, = 2.0(—) is 
too high now (though the sign is right) and (y2360)? = 2180 seems large. The 
4000 ev. point is now definitely too high. 
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The large reduced width of the 2360 ev. peak compared to the reduced widths 
of the first two resonances might lead one to suspect that the last peak is 
really two or more resonances which have not yet been resolved. However, till 
more experimental data are available on this point we can use the makeshift of 
adding a fixed AR,;*, representing high energy J = 3 resonances, to pull down 
the upper wing of the 2360 ev. peak. This will also raise the points of the low 
energy side, but not much. Since it is a fixed AR,,+, we must decrease (7337)? 
and (72360)? to get o,,(0) = 2.0, then adjust (y10s0)? to get agreement at 750 ev. 
The results are shown in Fig. 4, where we see that the over-all agreement is 
rather satisfactory, except that og, is again too large (2.0(—) instead of 
1.7(—)). 

The value of ARi,+ = 0.1610-'? cm. may well not be significant of itself. 
From the formulas [17] and [18] we see that o,, and ogo, depend only on the 
combination (a— Ri), so that adding a constant AR, to Ri: is equivalent to 
changing the radius a. Only the difference (a— AR,,) matters. From this point 
of view it is sensible to talk about two different equivalent radii a+ and a-, 
corresponding to different fixed AR,,+ and ARi,~. However, in this work it 
was not found necessary to use a fixed contribution from distant J = 2 
resonances. 

As a variant on Fig. 4, we try adding a J = 3 resonance somewhere above 
2360 ev.; this would pull down the 4000 ev. point and have much less effect at 
0 ev. than a constant AR,,+ does. There may be a resonance around 7000 ev. 
with (yz000)? = 1000. We raise (72360)? to 1900 to make o,,(0) = 2 barns, and 
decrease (71080)? to 66 to get a fit at 750 ev. We put a = 0.55 as usual, but for 
reasons to be explained later, Ty is raised to make on, = 14.3 at 2200 m./s. 
instead of 13.4. The result is Fig. 5, where the over-all agreement is seen to be 
excellent, except above 4000 ev. and for oon, which is 2.0(—) instead of 1.7(—). 

The reason that in Fig. 5, o,, was taken to be 14.3 at 2200 m./s. instead of 
13.4 can be seen from Fig. 6. This is a plot of the experimental low energy 
cross section against E~-?, where E runs between 1 and 0.02 ev. According to 
equation [15], this should be a straight line with y-intercept ¢,,(0) and slope 
proportional to Ty. However, only the lower part of this is a straight line; the 
upper part drops away. The straight line part has an intercept of 2.0 barns as 
it should, but is 0.9 barns above the experimental curve at 2200 m./s. Adjusting 
I, to make o,, = 13.4+0.9 = 14.3 barns at 2200 m./s. (instead of 13.4) hence 
gives a better fit in the energy region where the experimental curve of Fig. 6 is 
yet straight. 

The reason for the bow in the graph is not clear. A negative energy resonance 
(either J = 2 or J = 8) could cause such a fall off by decreasing Ri, and so 
decreasing o,,. However, the decrease would have to be of the order of 1 to 2 
barns, and since the total scattering cross section is only 2 barns, this seems 
improbable. The bow may be due to a crystal effect. In this case the curve 
should straighten out again as in the curves of Bendt and Ruderman (2). 

In considering negative energy resonances one must keep in mind that these 
decrease (g,Ri1t++g_Rit-), and if strong enough could make (a—g,Rut 
—g_R,-) positive at thermal energies, thus destroying agreement with the 
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Fic. 6. Low-energy total neutron cross section o; of Mn® as a function of E~+ (i.e, of 1/k, 
couueen: As in the other figures, the experimental points shown are intended only to indicate 


how much these scatter. According to theory (equations [15] and [16] of the text) the curve 
should be straight; actually only the lower part of it is. 


negative sign of on. For instance, using the parameters of Fig. 4, 
(a—g,Riut+—g_Rir-) is —0.4X10-? cm. at thermal energy, so any contem- 
plated negative energy resonance should not change g,Ri:1++g_Ri:- by more 
than 0.4X10-!? cm. 


CONCLUSION 


With no set of parameters were we able both to fit the cross section curve 
entirely and to account for the value of oo, accurately. However, rather close 
agreement was obtained with the parameters of Table I. Upon adding a fixed 


TABLE I 
RESONANCE PARAMETERS 
Ej, Yn’, Tn, 
ev. J ev-cm. ev. 
337 2 290 x 10-” 23 
1080 3 75X10-8 11 
2360 3 1850 x 10-” 395 





Tr, = 0.5 ev., a = 0.55 X10-" cm. 


AR,*+ = 0.16X10-!2 cm. to represent the contribution of distant, positive 
energy, J = 3 resonances, the experimental curve for the total cross section 
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and the scattering and absorption cross section at 0 ev. were fitted well. 
However, the value of the coherent scattering cross section at 0 ev. was not in 
agreement with experiment (calculated 2.0(—), observed 1.7(—) barns). In 
the thermal energy region, a not as yet understood effect causes the total cross 
section curve to deviate from the straight line one would expect. This may well 
be some crystal effect, in which case we may be misinterpreting the measured 
coherent cross section of 1.7 barns. 


Hence, the X matrix theory can fit the present experimental data quite well, 


but not perfectly. More comprehensive data on Mn®* are required before it 
can be stated definitely whether the discrepancies are real or not. 
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THE vz; INFRARED BAND OF AMMONIA! 


By C. CumMMING? 


ABSTRACT 


An absorption band in gaseous ammonia at 3 yu is observed at high resolution 
and attributed to the doubly degenerate fundamental v3. The constants vp = 3442.9, 
Ao = A; = 6.13, Aig = 0.29, and B; = 9.75 cm.— describe the band. Many of the 
observed lines are doublets. This doubling is a result of the inversion of the 
ammonia molecule and the value 1.1 cm.~! is calculated for the inversion 
doubling of the 0010 state. 


INTRODUCTION 


The fundamental band of the doubly degenerate »; vibration of gaseous 
ammonia has been observed in infrared absorption by Barker (1) and recently 
by Wood, Bell, and Nielsen (12). Neither investigation has led to an analysis 
of the band. A study of the Raman band has been reported (4) in which an 
analysis was carried out and the frequency of the band origin, 3444 cm.~!, 
was determined. Because of the comparatively low resolution of the Raman 
spectrum, due largely to pressure broadening, many of the measured lines were 
blends; accordingly the results of that analysis are less sure. It therefore 
seemed desirable to attempt an analysis of the infrared band. 


EXPERIMENTAL 


The band was observed in absorption using a Pfund-type spectrometer 
with an Anderson 600 lines/mm. grating and a lead sulphide detector. The 
practical resolution of the spectrometer at 3 » was about 0.3 cm.—! when the 
source of radiation was a Nernst glower. Because of pressure broadening in 
ammonia, pressures of less than 6 cm. Hg had to be used to take full advantage 
of the resolution of the spectrometer. 

Frequencies corrected to vacuum were found for the ammonia lines by 
comparison with known frequencies falling at the same grating angle. Certain 
lines in the 1.4 w band of water measured by Plyler, Gailar, and Wiggins (8) 
were used in the 3380-3500 cm.~! region, and the v3; band of hydrogen cyanide 
provided standards between 3300 and 3400 cm.—!. The rotational constants 
B’, B’’, D’, and D” for the latter band can be calculated very precisely from 
the results of Rank, Shearer, and Wiggins for 2v3; (9), and the v3 band origin 
3311.47 was fixed by comparison with Hg 10139.67 and Hg 15295.82. The 
accuracy of the ammonia frequencies is probably about 0.2 cm.—! and is 
limited by the interpolations between standards and by a slight nonuniformity 
in the rotation of the grating. 

Unfortunately atmospheric absorptions in the spectrometer obscured the 
spectrum almost entirely above 3500 cm.~'; almost half of the ammonia »; 
band could therefore not be measured. A number of weak lines below 3300 

1Manuscript received August 19, 1956. 
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cm.~! not attributable to v3 were observed. It is presumed that these lines belong 
to the 2v, band, but an analysis of these lines was not attempted. 


RESULTS AND DISCUSSION 


The frequencies of the measured lines are given in Table I. It will be noted 
that many of the lines are close doublets. This doubling is undoubtedly a 
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TABLE I 
S IN THE »3 INFRARED BAND OF AMMONIA 


Identification 


(J'K'-J"K") 


4,4-3,3 
2,1-1,0 
2,2-1,1 


2,2-3,3 
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3,2-4,3 
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Identification 
(J R= SE") 


4,4-5,5 
3,1-4,2 


4,1-5,2 
5,3-6,4 


6,5-7,6 


splitting due to inversion and will be discussed later. Since an insufficient 
number of the doublets have been resolved to justify analyses of the separate 
subbands, the unresolved lines and mean frequencies of doublets are used to 
calculate molecular constants which would describe the vibrational states in 
the absence of inversion doubling. 

The rotational levels for a symmetric top in an excited state of a degenerate 


vibration are, according to Teller (11), 
Fin (J, K) = Bey J (J+1) + (Aty— By) K? F 2ApafK, 


where ¢ is the vibrational angular momentum and the other symbols have their 
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usual significance. The vs band of ammonia corresponds to transitions to such 
a set of levels from the ground state rotational levels 


Fio(J, K) = BoJ(J+1)+ (Ao—Bo)K? 


in accordance with the selection rules AK = +1, AJ = 0, +1. Although 
there is the further restriction that only the +/ levels combine with the ground 
state for AK = +1, only the —/ levels for AK = —1, the band nevertheless 
has a fairly complex structure. To facilitate the identification of the observed 
spectrum, the intensities of the allowed transitions were calculated using the 
H6nl-London formulae (7) and the known populations of the ground state 
rotational levels. When the frequencies of the most probable transitions in the 
infrared band were calculated using the constants* from the Raman band (4), 
most of the lines in the observed spectrum were readily identified. 

The constants vo, A, Ao, Bi, and Ai¢ were recalculated making use of the new 
data. Regrettably the data are not extensive enough to allow a determination 
of the centrifugal stretching terms —D,J?(J+1)?—DjxK°J(J+1)—DgK‘4 
which would be expected to have an appreciable effect on the spectrum of such 
a light molecule. In carrying.out the analysis, therefore, corrections were 
applied to the measured frequencies to eliminate centrifugal stretching using 
the constants calculated by Slawsky and Dennison (10) as modified by 
Hansler and Oetjen (6). These are D; = 8.1X10-4, Dye = —15.5X10-4, 
De = 6.3X%10-* cm... 

Most of the strong lines fall into three series which can be designated 
PP(J,K = J), ’P(J, K = J—-1), and ’P(J, K = J—2). The frequencies of 
these lines can be represented by the equations 


PP(J, K = J) = [vo+A1—2A1¢—B,]—K[2A 1—2A 1 + Bo — Bi] +K2[A1—Aol, 

PP(J, K = J-1) = [v+Ai1—2A1¢—B,—2By]—K[2A1—2A ¢ +3Bo—3B;] 
+K?*[A,—Aol, 

PP(J, K = J—2) = [v+A1—2A+Bi1—6Bo]—K[2A1—24 ¢ + 5Bo—5Bi] 
+K*[A,—Apl. 


Since the graphs of these frequencies plotted against K form straight lines, one 
concludes that A»— A, is small. If Ag = Ai is assumed, then the intercepts and 
slopes of these graphs can be used to calculate certain molecular constants. 
The differences in the slopes give directly By —B, = 0.19, and since By = 9.94 
is well known (4), B; = 9.75. The intercepts and slopes allow one to calculate 
only the combinations v»»—Ai¢ and A; — Aig; in order to determine the 
constants vp, Ai, and ¢ independently, recourse was had to the frequency 
interval between the lines of the “S(J, K = J) series of the Raman band (3), 
which gave 


The constants vp = 3442.9, A; = 6.13 = Ao, and ¢ = 0.047 cm.—! were thus 
found. The spectrum calculated on the basis of these constants is shown below 
the observed spectrum in Fig. 1. 


*yo = 3444, B, = 9.78, Ao = Ai = 6.24, Aig = 0.28 cm. 
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Fic. 1. Tracing of the v; infrared band of ammonia. Superposed are the Q and R branches of 
v, (shaded area) and atmospheric absorptions (dotted line). The strongest lines in the vs band 
are plotted under the tracing, with frequencies calculated from the results of the analysis. 


The new values of the molecular constants are in substantial agreement 
with those determined solely from the Raman spectrum; this investigation 
thus confirms the interpretation of the Raman spectrum. In fact the differences 
between the new and earlier values can be attributed mainly to the new values 
for D; and D;x found by Hansler and Oetjen. The analysis can be criticized 
for its assumption of Ay—A:1 = 0, since this constant could be as large as 
+0.05 without manifesting a significant curvature in the plots. Evidently an 
accurate value for Ay and the moment of inertia of ammonia about the sym- 
metry axis must await more extensive data. 

In rotation—vibration transitions from the ground state of ammonia the 
symmetry selection rule for dipole transitions — <> +, and + «+> +, — «+!> — 
means that for parallel infrared bands with AK = 0 such as »; a doubling equal 
to the sum of the inversion splittings of the initial and final states should 
occur; on the other hand, for perpendicular bands with AK = +1 such as 7; 
the observed doubling should be the difference. The doubling of lines in the »; 
band was noted in the early investigation of Dennison and Hardy (5) and can 
also be seen in Fig. 1; apparently the smaller doubling in the »; band has 
hitherto escaped detection. 

The splitting in the 0010 state can be estimated from the data since the 
doubling in the ground state is known from the microwave spectrum (2). 
Fig. 2 is a graph of inversion splittings in the 0010 state for the rotational 
levels (J, K = J), plotted against J. By extrapolating to J = 0, one finds the 
value 1.1 cm.~! for the inversion splitting, which is not markedly different from 
the corresponding splitting (0.7935) for the ground state (2). Although the 
splitting due to inversion in ammonia has been evaluated theoretically only 
for the verve series of states, one would expect the splitting to have little depen- 
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Fic. 2. Determination of the inversion splitting for the 0010 state of ammonia. The graph 
shows the inversion splittings of the (J, K = J) levels plotted against J. 


dence on v3, since there is at most a second order change in the height of the 
NH; pyramid for this vibration. The large dependence of this doubling on 
J and K is not so readily understood. 
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THE FISSION YIELDS OF THE CESIUM ISOTOPES FORMED 
IN THE THERMAL NEUTRON FISSION OF U?® AND THE 
NEUTRON ABSORPTION CROSS SECTION OF Xe’! 


By J. A. PetrusKa,? E. A. MELAIKA, AND R. H. TOMLINSON 


ABSTRACT 


The absolute fission yields of the cesium isotopes occurring in the thermal 
neutron fission of U* have been determined with a mass spectrometer using 
isotope dilution techniques. Values of 6.59%, 6.41%, and 6.15% have been 
obtained for the yields of Cs!*3, Cs!%5, and Cs!8’ respectively. The neutron absorp- 
tion cross section of Xe!5 has also been measured to be 866 times greater than 
that of B!° for a Maxwellian distribution of neutron velocities corresponding 
to a neutron temperature of 57°C. This ratio gives a thermal neutron absorption 
cross section of 3.47 X 108 barns for Xe!*5 assuming it is a 1/v detector. 


INTRODUCTION 


Neutron absorption by the fission product isotope Xe'*> reduces the re- 
activity of nuclear reactors, particularly when they are operated at high 
thermal neutron fluxes. Reliable values for the cross section and the fission 
yield of this isotope are therefore of considerable interest in pile operations. 
Since the primary yield of Cs'* is probably small, the yield of Xe!** may be 
estimated from the cumulative yield of Cs!*5. The relative yields of the cesium 
isotopes in U**5 fission have been obtained from mass spectrometric data by 
Inghram, Hess, and Reynolds (8) and Wiles, Smith, Horsley, and Thode (14). 
In the former work the relative yield of the mass 135 isotope appears to be in 
considerable error owing to the failure of these authors to account adequately 
for the neutron capture in Xe!*, In the more recent work of Wiles et al. the 
relative yield for mass 135 was obtained by measuring the isotopic ratios of 
fission products formed under conditions where the neutron absorption by 
Xe!* was negligible. Since the large resonance absorption of neutrons by Xe!*5 
occurs in the thermal energy region, it was possible to minimize this absorption 
by shielding the uranium with cadmium during the irradiation. 

Absolute fission yields for each of the cesium isotopes may be obtained from 
these relative yields with the aid of a suitable normalization factor. Katcoff 
and Rubinson (9) have given a value of 6.62% for the yield of the mass 133 
fission chain. On the other hand the yield of this mass chain may be estimated 
to be 7.25% from the 3.2% yield for I'*! given by Bartholomew, Brown, 
Hawkings, Merritt, and Yaffe (1) and the 2.265 value for the ratio of mass 
133 to mass 131 given by Thode et al. (14). It is therefore of some interest to 
establish which of the values for the cumulative yield of the mass 133 chain 
is most reliable since this not only fixes the absolute values of the cesium 


1Manuscript received June 9, 1955. 
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3Vaffe and Brown have given the value of 3.1% for I'*! relative to a 6.1% yield of Ba'*°. The 
yield of Ba‘ has recently been given by Yaffe, Thode, Merritt, Hawkings, Brown, and Bartholomew 
(16) as 6.82% which therefore gives 3.21% for the yield of I'. A further revision of these values 
causing a 3-4% increase may also be required if a self-shielding correction is applied to the Ba'*® 
yield. 
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yields but also those of the xenon isotopes whose relative yields in fission 
products have been determined. The thermal neutron absorption cross section 
of Xe!*> has been reported to be 3.5 X 10° barns, but no details of this determina- 
tion have been released (10). 

In the present work the cumulative yield of Cs'** has been determined from 
measurements of the number of cesium atoms produced during a given irradia- 
tion and a knowledge of the extent of the U?* fission. The number of cesium 
atoms produced was measured with the mass spectrometer using isotope 
dilution techniques. The number of U** fissions was calculated from the change 
in the B!°/B" ratio in a simultaneously irradiated boron flux monitor and the 
known ratio of the neutron fission cross section of U*** to the neutron absorp- 
tion cross section of B!*. These data also give a direct determination of the 
yield of Cs'8’, but not Cs'*5, since the neutron flux for this irradiation was 
large enough that a considerable fraction of the 135 mass chain was lost 
through neutron absorption by Xe'**. To obtain the yield of Cs!** relative to 
Cs!83 a second sample of natural uranium was irradiated in a thermal column 
where the neutron flux was sufficiently low that neutron absorption by Xe!*® 
was small. 

The relative yields from these two samples also made possible an evaluation 
of the relative neutron absorption cross sections of Xe!*® and B’® since the 
decrease in the cumulative yield of the 135 mass chain owing to neutron 
absorption by Xe!*® may be related to the concomitant change in the B!°/B" 
ratio of the associated flux monitor. 


EXPERIMENTAL 


Two samples of natural uranium metal were irradiated with thermal 
neutrons in the N.R.X. (Chalk River) reactor. Sample A weighing 2.966 gm. 
(a rectangular parallelepiped of approximate dimensions 0.8 cm. X0.4 cm. X 
0.5 cm.) along with a quartz container filled with about 1 ml. of BF; gas at 0.5 
atm. pressure was irradiated with thermal neutrons for a 36-day period. 
The neutrons were considered to have a Maxwellian distribution of velocities 
corresponding to a temperature of 57°C. Sample B consisted of several similar 
samples weighing approximately 20 gm. each. These were simultaneously 
irradiated for a period of 80 days in the thermal column. 

In order to minimize contamination of the fission products by natural 
isotopes, polyethylene or teflon containers and pipettes were used in the 
preparation of reagents and the separation of the fission products. The 
uranium was dissolved in freshly distilled nitric acid and, where isotope dilu- 
tion was required, known amounts of natural Cs (Cs!**) were added to aliquots 
of this solution. The uranium was separated from the fission products by 
precipitation of uranium tetroxide from a neutral solution of the nitrate using 

- hydrogen peroxide as precipitant. The supernatant solution containing most 
of the fission products was concentrated by evaporation and transferred 
directly to the ion source filament of a 90° sector mass spectrometer. Separa- 
tion of the fission products from each other was unnecessary since ion emission 
of cesium occurs at a much lower temperature than any of the other elements 
which produce ions in this mass region. 
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The abundances of the boron isotopes in the BF; used for a flux monitor 
were determined with a 180° mass spectrometer. The ratios of the boron 
isotopes in the unirradiated BF; and the irradiated BF; were compared 
directly. Since both samples came from the same supply, the original isotopic 
composition of the irradiated BF; was identical to that of the unirradiated 
BF;. 


RESULTS AND CONCLUSIONS 


Thermal Neutron Fission Yields of the Cesium Isotopes 


The observed mass spectrometric ratios for the long-lived and stable isotopes 
of cesium occurring in the thermal fission of U?** as well as values of the ratios 
obtained after isotope dilution are shown in Table I. The errors shown for the 


TABLE I 


THE RATIOS OF CESIUM ISOTOPES OCCURRING IN THE THERMAL NEUTRON FISSION OF U2 
TOGETHER WITH ISOTOPE DILUTION DATA 


Isotope 


Cs133 Cs} Cs137 


Sample A 








.768+0 .003 








Mass spectrometric ratio* 1.074+0.003 0 
Ratio after isotope dilution 
I. 2.58410" atoms Cs!*8 added 
per gram of uranium 1.664 0.768 1.000 
IT. 3.117 X10" atoms Cs!*3 added 
per gram of uranium 1.789 0.768 1.000 
Atoms of fission Cs!3 per gram of 
uranium X10 
I 4.702 
I] 4.680 
Average 4.653+0.011f 
Sample B 
Mass spectrometric ratio* 1.072+0.002 1.034+0.002 1.000 


Ratio after isotope dilution 
I, 2.013 X10" atoms Cs!53 added 
per gram of uranium 1.937 1.034 1.000 
I]. 2.18410" atoms Cs!¥3 added 


per gram of uranium 2.003 1.034 1.000 
Atoms of fission Cs'3 per gram of 
uranium X10" 
I 2.488 
II 2.515 
Average 2.50 +0.01 


* Corrected for 26.6-year half-life of Cs'8* (18). 
{Corrected for 0.8% contribution of Cs! from fast fission of U*%, 


mass spectrometric ratios of Sample A are the standard deviations of the 
observed values from the average for a large number of spectrograms on this 
single sample. For Sample B the value given represents the average of 11 
determinations similar to that of Sample A made on fission products separated 
from three different samples of uranium. The precision shown for Sample B is 
the standard deviation of these 11 determinations from the average. Table I 
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also shows the number of Cs!** atoms produced per gram of uranium calculated 
from the fission product and isotope dilution ratios. The values obtained 
include Cs'** produced by the fast fission of U?*8. This contribution has been 
estimated to be 0.8% for Sample A and the final value shown in Table I has 
been corrected to give only the number of Cs!** atoms produced in the thermal 
neutron fission of U**5, 

From the number of cesium atoms produced per gram of uranium it is 
possible to determine the fission yield if the number of fissions per gram is 
known. This may be calculated if the thermal neutron flux in the uranium is 
known. The flux associated with Sample A was monitored with BF; which 
showed a ratio of B!°/B!! before irradiation to that after irradiation of 
1.029+0.001. Thus, if the only change in the boron isotope ratio was due to 
the loss of B® by the B!°(m, a)Li’ reaction, then 


1.029 = exp 50! °¢t, 
the neutron absorption cross section of B’°, 


where — 5a! 


? 
t 


the neutron flux in the boron, 


the time of irradiation. 


From this relation it may be shown that 
(1] soot = 0.0286. 


When the depletion of U*** is small, the total number of thermal fissions 
produced per gram of natural uranium is given by the expression 


[2] N;= 9207°°h'tNoas, 
where N, = the number of thermal fissions of U*** per gram of natural 
uranium, 
9204°° = the thermal neutron fission cross section of U**, 
¢’ = the thermal neutron flux in the uranium, 


No35 = the number of U?** atoms per gram of natural uranium. 


In order to evaluate equation [2] with the aid of equation [1], it is necessary 
to know the relationship between the flux in the boron monitor and the flux 
in the uranium. The flux reduction in the uranium or the self-shielding 
correction may be evaluated from the data of Case, Hoffmann, and Placzek 
(3). These authors have shown that, if an absorbing object is placed in a 
uniform isotropic neutron bath, the ratio of the average neutron flux within 
the object to the neutron flux at a distance outside the object is equal to Po 
which is a parameter known as the escape probability. Values of the collision 
probability, P. (which is equal to 1—Po), have been tabulated for spherical 
objects where the ratio of the object radius to the absorption mean free path 
is known. If the object placed in the isotropic neutron bath may scatter as 
well as absorb neutrons, the escape probability is reduced. Case, Hoffmann, 
and Placzek (3) have shown, however, that scattering must predominate to a 
very great extent before the escape probability is materially affected. For 
natural uranium, where the neutron absorption and scatter cross sections 
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are approximately equal (10), the effect of scatter may be neglected in the 
evaluation of the self-shielding correction. 

The absorption mean free path of the neutron may be calculated from the 
relation 
[3] h = 1/Noett, 


where \ = the absorption mean free path, 
N = the number of atoms per cm.’, 
Jett = the effective neutron absorption cross section. 


It has been assumed that the neutron absorption cross section of natural 
uranium varies inversely as the neutron velocity and that the neutron velocity 
distribution in the reactor was Maxwellian, and therefore the effective neutron 
absorption cross section of the uranium may be calculated from the relation 


Va, /293 
920ett = 920ab ~~ 5 4 - ’ 


where 920,» = neutron absorption cross section of natural uranium for 2200 
m./sec. (neutrons) 
= 7.42 barns (10), 
T = neutron temperature 
= 57°C. for experimental conditions. 


Thus the effective neutron absorption cross section of natural uranium is 

6.20 barns which when substituted in equation [3] gives: 
0.6023 X 10° XK 18.7 X 6.20 X 10 

That is, the absorption mean free path of neutrons having a Maxwellian 
distribution of velocities corresponding to a neutron temperature of 57°C. is 
3.42 cm. in natural uranium. If Sample A weighing 2.966 gm. were spherical 
in shape, the radius of this sphere would be 0.336 cm., and the ratio of this 
radius to the mean free path would be 0.0985. From the tables given by Case, 
Hoffmann, and Placzek (3), a value of 0.070 is obtained for P,,. Since the actual 
shape of the uranium in Sample A was a rectangular parallelepiped (approxi- 
mately 0.8 cm. X0.4 cm. X0.5 cm.), the flux reduction would not be as great 
as that calculated for a sphere. This factor, however, is not very critical unless 
the surface to volume ratio is much greater than that of a sphere. For example, 
Case et al. (3) have shown that a spheroid with major to minor axis ratio of 
5/3 would have a value of P, equal to 0.98 that of a sphere of equivalent 
volume. For Sample A therefore this shape factor gives P, = 0.068 and hence 
the self-shielding correction Py = 0.932. Calculations show that the accumula- 
tion of fission products such as Sm"° and Xe!'*® would introduce less than 
0.1% increase in this value for the particular irradiation conditions. Equivalent 
calculations show that the self-shielding correction for the boron flux monitor is 
negligible. The whole self-shielding correction for the uranium in relation to that 
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of the boron is therefore applicable providing that the flux in the boron was not 
reduced by the presence of the uranium. Since the relative positions of the 
uranium and the boron were not recorded, exact evaluation of this influence 
can not be calculated. There was, however, no scattering material close to or 
between the uranium and the boron and hence the flux reduction in the boron 
was merely the ‘‘shadow”’ cast by the uranium on the boron. It is estimated 
from the sizes of the uranium, boron, and the irradiation capsule that nat more 
than 10% of the neutrons passing through the boron had passed through the 
uranium. Since the average flux reduction in the uranium itself as calculated 
above was 6.8%, the reduction of the flux in the boron resulting from the 
presence of the uranium would have been less than 1%. The effective flux 
within the uranium is therefore taken as 0.932 that in the boron, i.e. ¢’ = 0.932 ¢ 
for Sample A. This relation, together with equations [1] and [2], gives 
235 


[4] N, = 0.9320.0286 x 224 Noss. 


50 





Yaffe et al. (16) have given a value of 
ono; /so° = 0.1456 
for identical neutron irradiation conditions and 


Nose = £:023X 10" 0.00719 
— 238.07 ; 


where the abundance of U?*5 in natural uranium is 0.719% (4). Substitution of 
these values in equation [4] gives NV; = 7.06X10'* thermal neutron fissions 
per gram of natural uranium in Sample A. From the number of Cs!** atoms 
produced per gram of uranium in Table I a value of 6.59% is obtained for the 
fission yield of Cs'%*. In view of the many parameters which have entered the 
calculation of this fission yield, it is difficult to assess the absolute accuracy 
of the value. The duplicate isotope dilutions and the mass spectrometry 
of the cesium isotopes would probably introduce less than 1% error. The 
largest apparent source of error is in the flux determinations from the change 
in the boron isotope ratio. The observed change was small and the precision 
of this value was only +3.5%. Other sources of error include the self-shielding 
correction, cross sections, etc., which could lead to a revision of the 6.59% 
yield if more reliable data and methods of calculation become available. 
The yield of Cs'*7 may be obtained from the absolute value of the Cs!*? 
yield and the relative yield of Cs!** to Cs!8? for Sample A. However, the 
extent of neutron absorption by Xe'*> was sufficiently large for this sample 
that it is difficult to evaluate the yield of the 135 mass chain. More reliable 
values of the relative yields of the cesium isotopes can be obtained from the 
data for Sample B. In this case the correction to the 135 mass ratio for neutron 
absorption is small compared to that for Sample A and the yield may therefore 
be estimated with greater accuracy. Wanless (13) has shown that the ratio 
of the observed yield of Cs!** to the cumulative yield of the 135 mass chain is 


given by equation [5]: 








646 CANADIAN JOURNAL OF PHYSICS. VOL, 33 


N 135 J» am Kt 1 don » soni 
Bl itm | a rg te YP 1 oR 


Y135 
_exp( ey 
~K(K—s3\"*) 
where N,;; = the observed number of Cs!*° atoms produced per gram of 
natural uranium, 
Yi35 = the cumulative yield of the 135 mass chain in atoms per gram 
of uranium, 
53\/%5 = the decay constant for I'** = 0.103 hr.— (6), 
54/85 = the decay constant for Xe!*5 = 0.0752 hr.—! (2), 
s47'*5 = the thermal neutron capture cross section of Xe!*, 
es 5Atd1354 gol 85Q’, 


Equation [5] may be approximately represented by equation [6] if the time of 
irradiation is long compared to the half-life of Xe**. 


Ni35 ao 

(6) Yiss sane og 

From equation [6] it may be seen that for a sufficiently low value of the flux 
the observed Cs'** yield is equivalent to the cumulative yield of the 135 mass 
chain. In the case of Sample B this condition is approached and therefore the 
ratio Ny35/Y135 may be evaluated with accuracy even if considerable error is 
involved in the values of the parameters required in the solution of equation 
(6). 

The neutron flux was not monitored for Sample B, but the isotope dilution 
data which give the number of cesium atoms produced may be used with the 
fission yield established from Sample A to evaluate the flux. Thus if the deple- 
tion of U?** is small, equation [7] gives the cumulative yield of Cs'**: 


[7] Nisz = Noss voy Ot fiss, 


where N33; = the number of Cs'** atoms produced per gram of uranium, 
fiss = the fractional fission yield of Cs'**. 


Equation [7] when rearranged and multiplied by 540'*® gives equation [8]: 


135 
ithe Ni33 540 

[8] so Gf = Nose fran * ooo’ 
23st fiss © 9207 


As an approximation the ratio s40'*°/9207°° may be considered independent of 
neutron velocity distribution and temperature and hence it may be evaluated 
using s4a!*5 = 3.5 10° barns and 920795 = = 549 barns (10). With the value of 
N33 for Sample B taken from Table J, f13; as determined for Sample A, and 

= 1947 hr. equation [8] can be solved: 

135 4+ 14 238.07 3. 5x 10° 
nob! = 2.50X 10" XE 993 x10" X0.00719X 1947 X0.0659* 549 
= 0.00068 hr. 
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Substitution of this value in equation [6] shows that the observed Cs'*® yield 
is 99.1% of the cumulative yield of the 135 mass chain for Sample B. It also 
may be seen that for the irradiation conditions of Sample B, a 10% error in 
estimating the product s4o1**’ will only introduce a 0.1% error in the cal- 
culated relative yield. The final value of the relative yield of Cs!*> shown in 
Table II has been obtained from the data for Sample B currected in this 


TABLE II 


THE RELATIVE AND ABSOLUTE YIELDS OF CESIUM ISOTOPES OCCURRING IN THE THERMAL NEUTRON 
FISSION OF U2% 


Absolute fission yield 


Isotope Relative ee 
yield This work Katcoff and Rubinson 

Cs133 1.072 6.59 6.62 

Cs185 1.043 6.41 6.40 

Csi87 1.000 6.15* 6.07* 


* Not corrected for delayed neutron emission. 


manner. The corrected relative yields in Table II have been normalized to the 
measured 6.59% absolute yield of Cs'** and shown as absolute thermal neutron 
fission yields. The yields obtained by Katcoff and Rubinson (9) based ona 
measured yield of Xe'** and the relative yields of Wiles, Smith, Horsley, and 
Thode are shown in Table II for comparison. The 6.59% yield of Cs!** is in 
excellent agreement with the 6.62% yield obtained by Katcoff and Rubinson 
(9) and suggests a possible error in the yield of 7.25% obtained from the I!*! 
yield and the ratio of the yields of mass 133 and mass 131. This latter ratio 
has recently been redetermined by Fleming and Thode (5) and the result is 
substantially the same as that previously published. The existence, however, 
of the 2.3-day Xe!**" isomer not known when the earlier ratio was determined 
necessitates a small correction to the Xe!**/Xe!'* ratio. It is therefore suggested 
that the I'*! yield may be about 2.9% instead of the currently accepted value 
of 3.2% obtained from the data of Yaffe et al. (1) when normalized to a 6.3% 
yield of Ba'*°. The difference in the mass 137 yields shown in Table II is 
largely a result of the different half-lives used to correct for the decay of Cs'%’. 
In the present work 26.6 years (15) has been used rather than the previously 
accepted value of 33 years (14). The absolute yield of Cs'*” shown in Table II 
is not the cumulative yield of all isotopes of mass 137 which occur in the 
thermal neutron fission of U?**. Delayed neutron emission by Xe!*? (with 22.0 
sec. half-life of the parent I'8”) results in a branching of the 137 mass chain. 
Pappas (11) has estimated the loss to the 137 mass chain to be 0.42% and 
hence the cumulative yield of isotopes of mass 137 is 6.57% providing no other 
chain branching processes occur. 


The Neutron Absorption Cross Section of Xe'** 


The value of the neutron absorption cross section of Xe!*® may be obtained 
from the measured decrease in the Cs!**/Cs'*7 ratio for Sample A compared to 
the fission yield ratio and the change in the boron isotope ratio of the flux 








648 CANADIAN JOURNAL OF PHYSICS. VOL. 33 


monitor. The value obtained is therefore independent of the measured 
absolute fission yields. With equation [5] it is possible to calculate a series of 
Ni35/Yis5 ratios for various values of the product 540!*°¢’. The value of 
540'*°o’ which corresponds to the observed Nj35/V135 ratio may be used in the 
evaluation of 540%, 

Unfortunately the 36-day irradiation of Sample A was not continuous. The 
reactor operated at a constant power output for a total of 650 hr. which in 
fact represented six separate irradiation periods of 99.8 hr., 111.6 hr., 104.0 hr., 
123.8 hr., 139.3 hr., and 71.5 hr. respectively. The periods of irradiation were 
separated by shut-down periods of 56.0 hr., 48.5 hr., 35.5 hr., 26.0 hr., and 25.5 
hr. respectively. The periods between irradiation, however, were sufficiently 
long compared to the half-life of Xe1** that it was possible to consider that each 
of the six irradiation periods was independent. It was therefore necessary to 
obtain six calculated values of Ni35/VY135 from equation [5] for an arbitrarily 
chosen value of 5401*°¢’ corresponding to each of the six irradiation times. 
These ratios were then weighted in accordance with the fraction of the total 
irradiation period which they represented and averaged io give a value of 
Ni35/ Y135 corresponding to the chosen value of 540!*°¢’. Graphical representa- 
tion of Ni35/ Y135 for various values of 540!*°g’ made possible the evaluation of 
547/*°’ which corresponded to the experimental 4435/ Y135 ratio. The experi- 
mental value of Ni35/ Yiss for Sample A may be obtained from the product of 
the observed Cs!*5/Cs!87 ratio (Table I) and the fission yield Cs'*7/Cs!'** ratio 
shown in Table II, i.e. Ni35/Yiss = 0.768/1.043 = 0.736. This value corre- 
sponds to 5401*°¢’ = 0.0355 hr.-! which with equation [1] and the relation 
¢’ = 0.932¢ gives 5401*5/50!° = 866+40 for the assumed Maxwellian distribu- 
tion of neutron velocities corresponding to a neutron temperature of 57°C. 
This ratio gives an effective neutron absorption cross section of 2.90 10° 
barns for these experimental conditions if the absorption cross section of 
natural boron is 755 barns (12) and the abundance of B!° is 18.83% (7). 
If it is assumed that the neutron absorption cross section of Xe'** varies 
inversely as the neutron velocity, the neutron absorption cross section of 
X: ° for 2200 m./sec. neutrons becomes 3.47 X10® barns, which is in close 
agreement with the value of 3.5 10° barns which was released jointly by the 
United States, the United Kingdom, and Canada (10). 
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GAUGE-INVARIANT FORMULATION OF 
QUANTUM ELECTRODYNAMICS! 


By P. A. M. Dirac 


ABSTRACT 


Electrodynamics is formulated so as to be manifestly invariant under general 
gauge transformations, through being built up entirely in terms of gauge-invariant 
dynamical variables. The quantization of the theory can be carried out by the 
usual rules and meets with the usual difficulties. 

It is found that the gauge-invariant operation of creation of an electron 
involves the simultaneous creation of an electron and of the Coulomb field 
around it. The requirement of manifest gauge invariance prevents one from using 
the concept of an electron separated from its Coulomb field. 


GAUGE TRANSFORMATIONS 


A gauge transformation of the electromagnetic field changes the potentials 
A*(u = 0, 1, 2, 3) according to the law 


[1] A* — A+ dS/dx,. 


In the usual formulation of electrodynamics, one restricts the potentials by 
the condition 
[2] 0A,/dx, = 0. 


Then, in order that the transformation [1] shall preserve this condition, it is 
necessary that 
[3] 0°S/dx, dx" = 0, 


so that the gauge transformations are restricted. 

It is possible to do without the condition [2], in which case S can be an 
arbitrary function of the four x’s. One then gets an electrodynamics with 
general gauge transformations. Such a formulation of electrodynamics has 
some advantages in that (i) it is simpler, from a certain point of view, as it 
does not make use of an unnecessary condition, and (ii) it imposes stricter 
limitations on any alterations one may make when one tries to improve the 
theory. 

It will be shown here how one can develop electrodynamics without the con- 
ditions [2], [3] and can quantize it, keeping all the time to a formulation for 
which the invariance under general gauge transformations is manifest. 


THE LAGRANGIAN AND HAMILTONIAN 


We must work with the dynamical variables on a space-like three-dimen- 
sional surface in space-time and see how they are connected with the dynamical 
variables on a neighboring surface. These connections constitute the equations 
of motion. For simplicity we shall here deal only with the flat surfaces x») = 

1Manuscript received August 23, 1956. 

Contribution from St. John’s College, Cambridge, England. 


This paper formed the subject of two lectures given at the National Research Council at Ottawa 
on 10 and 12 August, 1956. 
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constant, corresponding to instants of time. The generalization to arbitrary 
space-like surfaces can be made without difficulty, but requires elaborate 
notation. 

With the Lagrangian formulation, the variables that we work with are the 
coordinates and velocities of the field, namely A®, A’ (r = 1, 2,3), 0A°/Axo, 
0A‘/dxo, together with suitable variables to describe the electrons. The 
variables are to be taken at all points of three-dimensional space for a fixed 
xo, so they are functions of three parameters +x,. 

Let us see how the field variables are affected by a gauge transformation. 
We have 

A® + A°+dS/Adxo, A’ A'+0S/dx;, 


[4] | 0A°/dx9 > 0A°/dxo+0°S/dx0", 9A"/Ax9—> JA" /dxo+0°S/ Ax, Axo. 


If we take S and 0S/dx» to be zero at all points of three-dimensional space at 
the time considered, then there is no change in A°, A’, and 0A’/dxo. However, 
0°S/dxo? remains at our disposal, and so we can make an arbitrary change in 
0A°/dx9. This marks an important difference from the usual formulation, in 
which [3] holds, so that if S and 0.8/dxo vanish then 02S/0x 9? must also vanish. 
With the present theory, the way in which A° varies with the time is completely 
arbitrary, instead of being fixed by the equations of motion and initial conditions. 

The passage from the Lagrangian to the Hamiltonian can be made by 
standard methods. The work has been done previously (1) and need not be 
recapitulated, but the main features in it will be mentioned here to show up 
the differences from the usual theory. 

For electrons with spin, described by the four-component wave function y, 
interacting with the electromagnetic field, the Lagrangian is 


[5] L=J{RNih&PY+a, y’)+e P(A+a, A)y+m J om ¥—% Fur For} dex, 


where § means the real part of the term that follows, and the suffix 0 or r 
added to a field variable such as y means its derivative with respect to xo or 
x,. In the usual theory the last term of the integrand is replaced by 


(6] —4A,, A», 


The difference of the two terms, integrated over the whole of space-time, 
vanishes with the help of [2], and so the two Lagrangians give the same equa- 
tions of motion. 

The momenta B,(x; x2x3) conjugate to the coordinates A*“(x1 x23) are 
introduced in the usual way, as the coefficients of 5A*° in the variation of the 
Lagrangian. Now the Lagrangian [5] does not involve A”, and hence it makes 
By vanish. The vanishing of By does not produce a difficulty, but merely 
appears as a supplementary condition in the Hamiltonian scheme of equations. 
The other components B, of B are equal to F,», the components of the electric 
field. 

In the usual theory with the term [6] in the Lagrangian, By does not vanish 
and the four B’s are B, = —A,°. These B’s are not gauge-invariant quantities, 


as in the new theory. 
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The Hamiltonian that follows from the Lagrangian [5] is 
[7] H=f{-Rihkpa,p'—epa, A’y—mPany+h Fs F*— 4B, B}d*x. 
It is to be used in conjunction with the supplementary conditions 
[8] By = 0, 
[9] B,’+eWy = 0. 


The second of these is a consequence of the field equations, and since it involves 
only dynamical coordinates and momenta, it has to be taken into account in 
the Hamiltonian scheme. The supplementary conditions are all consistent, 
because their left-hand sides have zero Poisson bracket with one another and 
with the Hamiltonian. 

In the usual theory the last two terms of the integrand in [7] are replaced 
by 
[10] 3A ys A*'— 4B, Be 


and a further term must be added to the integrand, 
[11] —eAW, 


representing the potential energy of the electrons. The supplementary condi- 
tion [8] must be replaced by 
[12] Bo- , = 0 


which follows from [2], and the supplementary condition [9] still holds. 

An important feature of the new theory is that the variables A°, By do not 
appear in the Hamiltonian [7]. It follows that these variables are of no physical 
importance and can be dropped out from the theory. This is a simplification 
which does not occur in the usual theory, where A® and Bo appear in the terms 
[11] and [10]. 

The supplementary conditions can be looked upon as Hamiltonians giving 
rise to permissible changes in the dynamical variables at one instant of time. 
In the present case these changes are the gauge transformations. 

The supplementary condition [8] merely gives rise to an arbitrary change in 
A° with the other dynamical variables unchanged. This is a gauge transforma- 
tion [4] with S zero at the time considered, but not 0S/dxo. Such gauge trans- 
formations are automatically allowed in the theory when the variables A°, By 
are dropped out. 

To see the changes produced by the supplementary condition [9], let us 
work with the general linear combination of [9] at all points of three-dimen- 
sional space, say 


[13] G = fg(B,’+epp)d*x = 0, 


g being a general function of the three x’s. The supplementary condition [13] 
will make the dynamical variables vary with some parameter 7 according to 


dA'/dr = [A',G] = fg’[A’, B,* +e)'p']d*x’ 
[14] = fg'5(x—x’)"dix' = —g". 
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With the suffix a denoting one of the four components of y, 
dyq/dr = [a, G] = fe’ [va Bs’ +edo'yo’ |d*x’ 

—teh day fg’s (x—x’)Wp’ d3x’ 

—teh gq. 


¥ will vary with 7 according to the conjugate complex equation, and the other 
dynamical variables will remain unchanged. 

With 7 going from 0 to 1 say, [14] is just a gauge transformation [4], with g 
equal to —S at the time considered. Equation [15] shows how y changes with 
this gauge transformation. 


[15] 


THE GAUGE-INVARIANT QUANTITIES 

Our problem now is to develop electrodynamics so as to be manifestly 
invariant under gauge transformations of the type [14], [15]. We can accom- 
plish this by working entirely with gauge-invariant quantities. 

A gauge-invariant quantity satisfies the condition that it has zero Poisson 
bracket with the left-hand side of [13], or of [9] at all points in three- 
dimensional space. The electric field B, and the magnetic field F’* at all 
points are examples of such quantities. Besides these, there is a gauge- 
invariant quantity involving y, which we must now calculate. 

Consider the quantity 
[16] Wa" (x) = Ya(x)er™, 


where C is a linear function of the A’ over the whole of three-dimensional 
space, besides being a function of the three x’s. Thus C is of the form 
[17] C(x) = fc, (x, x’) A(x’) d*x’, 


with c,(x, x’) a number. By a suitable choice of c,(x x’) we can make the 
quantity [16] gauge invariant. 
We have 
[va*, G] = [War Ge’ + ips e' fc, (x x’)[A"(x’), G] d*x’. 


Using the results [15], [14], this becomes 
[y.*, G] = —ip,* {eh g+ Ser(x x')g’" d*x'} 
— ipa* {eh1g— fc, (x x')" g'd*x’}. 


Thus y,* is gauge invariant provided c,(x x’) satisfies the condition 

[18] Cr(x x’)? = eh 5(x—x’). 

This means that c,(x x’), considered as a function of the x’’s, is a vector whose 
divergence vanishes everywhere except at the point x’ = x, where it has a 


singularity of the form of a 6 function. 
Equation [18] has many solutions. One of them is the spherically symmetrical 


solution 


et) ae HO Me Xe 
[19] CGS) re ee 
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Any other solution differs from this one by a vector whose divergence with 
respect to x’ vanishes everywhere. Such a vector is of the form e,,,b'"’, 
where b’‘ = b'(x’) is some vector function of the x’’s, and e,;, has the value 
+1 for r,s, ¢ an even or odd permutation of 1, 2,3 and zero otherwise. With 
c,(x x’) changed in this way, the change in C is 


[20] Serer bd! Ade! = —f eye b!'A!* die’, 


which involves only the gauge-invariant quantities A"*—A*’ = F*’. Thus the 
various choices of c,(x x’) give various y*’s differing from one another by 
gauge-invariant factors involving the magnetic field. 

The variables y* with some choice of c,(x x’), together with their conjugate 
complexes * and the electric and magnetic field quantities, are the only 
independent gauge-invariant quantities. A gauge-invariant theory should 
involve only these dynamical variables and functions of them. This means that 
whenever y occurs, it should be accompanied by a factor e'° which makes it into 
y*. 

An important gauge-invariant quantity connected with the differentiation 
of Wis 
[21] (thy'+eA pet. 

To check its gauge invariance, we note that it is equal to 
ihy**+y* (hC'+eA*) = ihp**+y*f {hic,(x x’)' +0586 (x—x')} A’? dix’. 


The coefficient of A’’ here is a vector, with respect to the suffix 7, whose 
divergence with respect to the variables x’ vanishes, on account of [18]. 
Thus the integral is like [20] and is gauge invariant. It follows that [21] is 
gauge invariant. 

The gauge invariance of the Hamiltonian [7] now becomes manifest, as it 
may be written 
[22] H = f{—R*a,(iip’+eA’ pe!°—mp*an Ye +i FF + 4B, B,} dx, 
involving only quantities like [21] and other gauge-invariants. 

QUANTIZATION: GENERAL THEORY 

To quantize a dynamical theory, one must make the dynamical variables 
into linear operators satisfying the appropriate commutation or anticom- 
mutation relations, corresponding to their Poisson bracket relations. These 
linear operators operate on ket vectors, which represent the various physical 
states. 

A convenient way to describe the kets is to refer them all to one standard ket, 


|S) say, representing some standard state. The ket |S) is defined mathe- 
matically by certain linear conditions, say the conditions 


[23] mS) = 0, 


the parameter \ being introduced to distinguish the various conditions. The 
conditions [23] must be consistent, which requires that the commutator or 
anticommutator of two x«’s shall be a linear function of the x’s of the form 


[24] fos Ky dd 
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with the coefficients 2, on the left. Also the conditions [23] must be sufficient 
to fix |S) completely, except for a numerical factor. 

When the standard ket |.S) has been chosen, one can express any other ket 
|P) as the result of some linear operator ¢ applied to |S), 


[25] |P) = ¢|S). 


One may restrict ¢ to be a function of a certain complete set of commuting or 
anticommuting dynamical variables, in which case ¢ is called a wave function. 
However it is not necessary to restrict ¢ in this way. In the general case the @ 
corresponding to a certain |P) is not uniquely determined, as any two ¢’s 
differing by a quantity of the form [24] correspond to the same |P). 

If the classical theory contains some supplementary conditions 6 = 0, 
these appear in the quantum theory as restrictions on the kets |P) that 
represent physical states, 


[26] B|P) = 0. 
The standard ket |S) must satisfy these restrictions, 
[27] |S) = 0. 


This means that the operators 8 must be included in the set of operators xy. 
The restriction [26] becomes, with the help of [25], 


84|S) = 0 
or, from [27], 


(8, ]|S) = 0. 


Thus if we restrict ¢ to involve only quantities that commute with the §’s, 
the ket |S) necessarily satisfies the conditions [26] and so represents a physical 
state. It is sufficient to consider only such ¢’s to get all physical states. 


QUANTIZATION: ELECTRODYNAMICS 


In electrodynamics the natural way of choosing the standard ket |S) is to 
take it to correspond to the vacuum state, for which there are no electrons, 
positrons, or photons present. A difficulty arises here in that one does not know 
precisely what mathematical conditions [23] to take for this ket. 

The usual procedure is to make a Fourier resolution of the various field 
quantities y, Y, B,, F’* and interpret the various coefficients as emission and 
absorption operators of electrons, positrons, and photons. One then defines 
the vacuum state as the state to which the application of any absorption 
operator gives zero. The objection to this procedure is that the vacuum state 
defined in this way does not satisfy the equations of motion. If one starts off 
with this state, after a lapse of time there will be some electrons, positrons, 
and photons existing, forming what are usually called the vacuum fluctuations 
of the electromagnetic field. Physically this is not correct. The vacuum state 
ought to remain unchanged. 

The difficulty is the same in the usual formulation of electrodynamics as in 
the present gauge-invariant formulation. It must mean that the Hamiltonian 
is wrong in some fundamental way. No better Hamiltonian is known at 
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present. So we have the situation that we cannot define accurately the vacuum 
state. We therefore have to work with a standard ket |S) which is ill-defined. 
One can, however, do many calculations without using the accurate conditions 
[23] and the successes of quantum electrodynamics are obtained in this way. 

Let us ignore the difficulty and suppose we have a standard ket |S) repre- 
senting the vacuum state. Any other ket can then be expressed as |S), ¢ 
being some operator. We now take over the Hamiltonian [7] or [22] and the 
supplementary conditions [9] into the quantum theory. For this purpose 
we must express all the field quantities in [7] or [22] and [9] in terms of emission 
and absorption operators and then arrange that in any product containing 
both emission and absorption operators as factors, the emission operators are 
all to the left of the absorption operators. This rearrangement causes the 
Hamiltonian [7] or [22] and the left-hand side of [9] to change by certain 
infinite quantities, which will be denoted by ky and k(x). k(x) must be a 
constant while ky may also contain terms linear in the A’s. 

The supplementary conditions now read 


[28] (B,’+ely+k)¢|S) = 0 
for all x. They must be satisfied by the standard ket |S), so 
[29] (B,’+ef~+k)|S) = 0. 





We now see that [28] is satisfied provided ¢ commutes with B,’+eyy, that is, 
provided ¢ is gauge invariant. 

If @ involves only the electromagnetic field variables B,, F’*, the state ¢|.S) 
has only photons and no electrons or positrons. To get a state with one electron, 
we must take ¢ to be linear in the gauge-invariant operators y,*(x), thus 
[30] |S) = fza(x) pa* (x) dx |S), 
where the coefficients z,(x) may involve B,(x’), F’*(x’) for all values of x’. 
Similarly, to get a state with two electrons we must take ¢ to be quadratic in 
the y*’s, and so on. For positrons one must replace y* by #*. 

The representation of states in this way is not unique, i.e., corresponding toa 
given state there are many ¢’s. Without having |S) well-defined one cannot 
formulate precisely which ¢’s correspond to the same state. 


THE ELECTRON EMISSION OPERATOR 
The ket 
[31] ¥a*(x)|S) = Pa(x)e*|S) 


represents a state for which there is one electron, situated at the point x. 
Let us examine the electric field around this electron. We have 


BY Wa" |S)= Wa" Bi, |S)+ih[ Bi, va" ]|S) 
= y." Bi, |S)—hy,*[B’. C]|S) 
va" {Bi +h c,(x x") }|S) 





[32] 


from [17]. The physical interpretation of this result is that the value of the 
electric field B,(x’) for the state y,*|.S) exceeds its value for the vacuum state 
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|S) by hc,(x x’). If we take the spherically symmetrical solution [19] for 
c,(x x’), the excess electric field is just the Coulomb field of the electron at the 
point x, with the system of units here used, (the charge on the electron being 
+e). With any other solution for c,(x x’), the excess electric field is the Coul- 
omb field plus a field of pure electromagnetic radiation. ; 

We can now see that the operator ¥*(x) is the operator of creation of an 
electron together with its Coulomb field, or possibly the operator of absorption 
of a positron together with its Coulomb field. It is to be contrasted with the 
operator ¥(x), which gives the creation or absorption of a bare particle. A 
theory that works entirely with gauge-invariant operators has its electrons and 
positrons always accompanied by Coulomb fields around them, which is very 
reasonable from the physical point of view. 

By taking various solutions of [18] for c,(x x’) one can get the ket [31] to 
represent the state of one electron at the point x with its Coulomb field and 
with various possible distributions of photons. However, not all distributions of 
photons accompanying an electron can be described in this way, but only a 
certain simple class of them. To get a general photon distribution one must 
multiply [31] by a function of B,’ and F’” at all points x’. 

One might decide to work always with the spherically symmetrical c,(x x’) 
given by [19]. The convenient way of proceeding is then to resolve the electric 
field B, into a part consisting only of transverse waves, B', whose divergence 
vanishes, and a part consisting only of longitudinal waves, B!, whose curl 
vanishes. With the spherically symmetrical c,(x x’) it is easily verified that B‘ 
commutes with all the ¥*’s. Thus equation [32] becomes, with x and x’ inter- 
changed, 

pponoaestt 
3 


[33] B; va" |S) = va" Bits jsf |S). 


| 
x | 


If we could define the vacuum state so as to have no particles, we should 
expect it to have no longitudinal electric field, so that 


[34] B} |S) = 0. 


Then [33] reduces to 





[35] B! y* |S) = 72-75 yo |S) 


4m |x—x [> ** 


This means that the state y,’*|S), having an electron located at the point x’, 

is an eigenstate of the longitudinal part of the electric field at all points, the 

eigenvalues being just the Coulomb field. The result can easily be extended to 

any distribution of electrons. B! is always equal to the total Coulomb field. 
The electric energy term in the Hamiltonian becomes 


[36] 3 {B, B,d’x = 4 J (B; B,+B; B;) d’x, 
the cross terms vanishing. When this energy operator is applied to a ket 
¢|S) representing any distribution of electrons, with no positrons, the above 


result shows that the longitudinal part 4 /B? B' d*x is just equal to the energy 
of the Coulomb field of the electrons. This energy may be expressed, in inte- 
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grated form, as the sum of an interaction energy e?/4ar for each pair of elec- 
trons, together with an infinite self-energy for each electron. The theory now 
becomes identical with the usual quantum electrodynamics after elimination of 
the longitudinal electromagnetic waves. 

We can now see the significance of the usual method of elimination of the 
longitudinal waves. It implies that one looks upon each electron as always 
surrounded by a spherically symmetrical Coulomb field, and any departures 
from this field are ascribed to the existence of photons. The method provides 
us with a definite basis for making calculations, but one can see that it is not 
always suitable physically. 

Suppose for example that we are dealing with a state having an electron 
moving with a definite velocity in the absence of photons. (Velocity here 
means average velocity, without the ‘‘Zitterbewegung’’.) It would be reason- 
able to associate the electron with, nota spherically symmetrical Coulomb field, 
but a Coulomb field subjected to the Lorentz—Fitzgerald contraction appro- 
priate to its velocity. We should thus take for c,(x x’) the solution [19] dis- 
torted by the Lorentz—Fitzgerald contraction. With ¥* so defined and with 
the coefficients z,(x) in [30] chosen to give plane waves, we should have the 
best representation of the one-electron state that is possible with the present 
imperfect theory. It is an imperfect theory because the ket representing the 
one-electron state, like that representing the vacuum state, does not satisfy 
the equations of motion. 


THE MOMENTUM 
The total momentum M’ has the property that, for any field quantity 
E(x), 


[37] [é, M7] = &”. 
It is easily verified that if we take 
[38] M’ = [(B, A’ +Rihdy’)d'x, 


it has the desired property, for § = A‘, B,, y, or y, e.g., 
[A', M’] = J[A', Bi A +R in yp" J d's’ 
= [66(x—x’) A’ d’x’ = A”. 
This property is sufficient to fix M’ except for an arbitrary additive constant. 
The constant is fixed by the requirement that M’ applied to the vacuum state 


shall be zero. M’ is then completely determined. 
We must check that the expression [38] for M’ is gauge invariant. It may 


be written 


M’ = [(B, F°+B, A"+Rihpyp’)d'x 
[39] = f{B, FU+R* (ihy’+eA Yet }d'x—fA'(B,'+epp)d'x. 
The first of these terms is evidently gauge invariant. The second vanishes with 


the classical supplementary condition [9]. However, with the extra infinite 
constant term in the quantum supplementary condition [28] it does not vanish 
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and a difficulty arises. Evidently some change must be made in the theory, 
invoiving a more precise handling of the infinities, so that the second term 


vanishes also in the quantum theory. 

The first term in [39] then gives the total momentum. The part f B, F"* d*x 
just gives the momentum of the electromagnetic field, described by the 
Poynting vector. The remainder must be the momentum of the electrons. 


DISCRETE FARADAY LINES OF FORCE 
The equation [18] has a certain discrete type of solution which is of interest. 
Take 


feat x’) = —eh"e(x1—x1') 5(x2—x0") 5(x3—x3'), 
[40] [cole x’) = ¢3(x x’) = 0, 
where the function €(x) is 
e(x) = 1 forx > 0 
= 0 for x < 0. 


[40] is a solution of [18] and it gives for C a line integral 


[41] C= — ei f A’ (x; X2 %3) dx’. 


More generally, we may take for C any line integral 


[42] C= -aa” f A''dx,! 


_ 


extending from infinity to the electron point x. Such a C will make ¥*(x) = 
¥(x)e'© gauge invariant. 

Each choice of C corresponds to a choice of electric field around the electron. 
The choice [42] corresponds to the electric field being all concentrated in a 
single Faraday line of force, extending from the electron to infinity. The 
choice [41], which is more specific, makes the Faraday line of force extend in a 
straight line in the direction of the negative x-axis. 

One of the most fundamental unexplained features of present-day physics 
is the fact that electric charge always occurs in multiples of the electronic 
charge e. A reasonable attempt at an explanation could be founded on the 
assumption that quantization of the electromagnetic field results in the electric 
field being composed of discrete Faraday lines of force, each associated with 
the charge e, so that there is just one line ending on each charge e. The lines 
would then be the elementary concept in terms of which the whole theory of 
electrons and the electromagnetic field would have to be built up. Closed lines 
would be interpreted as photons and open lines would have their ends inter- 
preted as electrons or positrons. 

With such a theory it would be necessary, of course, to depart from the usual 
expression for the energy of the electric field, given by the last term in the 
integrand of [7], since this expression would have an infinite value for a 
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discrete Faraday line of force. We have seen that the Hamiltonian needs 
modification in any case, and we have here a possible basis for seeking a new 
one. A representation of the discrete lines of force could be provided by the use 
of operators y*(x) defined with C of the form [42]. 
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THE FLUCTUATING FIELD MODEL OF FERROMAGNETISM 
WITH PARTICULAR REFERENCE TO NICKEL' 


By F. D. Stacey 


ABSTRACT 


Following Néel’s calculations of fluctuations in intermolecular field, a model 
is suggested in which there is a mutual repulsion of elementary magnets, 
so that spontaneous magnetization is affected by a process of ordering of the 
magnets on the lattice sites as well as the ordering in orientation. Specific heat 
considerations indicate that, if it occurs, the second disordering process is also 
complete at the Curie point. To explain the observed spontaneous magnetization 
of nickel it is necessary to assume coupled electron spins and appreciable next- 
nearest neighbor exchange interactions. 


1. INTRODUCTION 


Néel (6, 7, 8) showed substantial correlation between the 
specific heat above the Curie temperature, 6,, of a ferromagnetic, the difference 
between this and the paramagnetic Curie temperature @,, and the curvature 
of the (1/x-7) graph above 6,, in terms of fluctuations of the intermolecular 
field, for which the same parameter was required to explain all the obser- 
vations. The calculations were based on two alternative models, neither of 
which bears close resemblance to the accepted structure of a ferromagnetic, 
but the meaning is clearly that elementary magnets being positive 3d holes 
are considered to have freedom of movement through a sea of 4s (conduction) 
electrons, but at any instant of time to have distinct positions. This approach 
has been criticized (12), being a contradiction of band theory in general and 
the collective electron theory of ferromagnetism in particular. Three develop- 
ments of Néel’s theory are needed: 

(i) Postulation of a more acceptable model in terms of recognized lattice 
structure. 

(ii) Introduction of band theory of which the success is too great for a 
serious model to ignore. 

(iii) Extension of the theory to explain the marked divergence at low tem- 

peratures of present experimental and theoretical spontaneous mag- 
netization curves. 


‘ 


‘anomalous”’ 


The present paper is an attempt to deal with requirements (i) and (iii) 
simultaneously, so that some conclusion regarding (ii) can be drawn from its 
success or failure. 

The model considered is a face-centered cubic lattice, of which only a known 
proportion of the lattice sites are occupied by elementary magnets. The pro- 
portion is 0.604 or 0.302 for the cases of the elementary magnets being single 
or coupled electron spins respectively. The magnets are allowed freedom of 
spontaneous transfer between lattice sites subject to the conditions that there 
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is no double occupation of any lattice site, and that the magnets experience a 
spin-independent mutual repulsion in addition to exchange interaction. This 
repulsion may arise electrostatically by virtue of the positive charges asso- 
ciated with the magnets. It is represented by mutual potential energies E, 
and £2 due to repulsions between pairs of magnets on nearest and next nearest 
neighboring sites. 

A thermal rearrangement of elementary magnets on the lattice sites causes 
changes in the numbers of nearest neighbors to the magnets and hence in the 
intermolecular field. The ordered and disordered limits of this thermal re- 
arrangement are considered, to determine the conditions under which the 
model is sufficient to explain the observed variation in intermolecular field in 
nickel. 

No further consideration is given to Néel’s problem of fluctuations, but the 
model is designed to give fluctuations in intermolecular field by virtue of 
fluctuating numbers of nearest neighbors. 


2. SPONTANEOUS MAGNETIZATION AND INTERMOLECULAR 
FIELD AS FUNCTIONS OF TEMPERATURE 

The quantum-corrected Weiss theory gives values of reduced spontaneous 
magnetization, ¢/o, as a function of reduced temperature, 7/6,;, which are in 
reasonable agreement with experiment, being closest for uncoupled spins 
(s = 3). Most tables of comparative values, such as Mott and Jones (4, p. 221), 
show that the disagreement is greatest at high temperatures, being about 10% 
at 7/0; = 0.9. Table I is designed to show that the major disagreement does 
in fact lie at low temperatures where small differences in o/go are critical. It 
is particularly important to make this comparison since Néel (8, p. 131) 
showed that the high temperature discrepancy disappears if 0, is used in 7/8. 
Table I gives values of intermolecular field which are needed to explain the 


TABLE I 
TEMPERATURE DEPENDENCE OF INTERMOLECULAR FIELD IN NICKEL 


Lr £ 6) (F) 
6 go Fo s—} Fo s=1 





0 1.0 0.115* 
0.05 0.99875 0.185 0.246, 
0.1 0.9960 0.312 0.369 
0.2 0.9875 0.513 0.593; 
0.3 0.9766 0.682 0.7734 
0.4 0.9585 0.807 0.895 
0.5 0.934 0.904 0.989 
0.6 0.892 0.963 1.033 
0.7 0.830 1.002 1.055 
0.8 0.737 1.025 1.061 
0.9 0.582 1.030 1.047 
0.95 0.451 1.024 1.034 
1.0 0.0 1.0 1.0 





*Extrapolated graphically. 


observed spontaneous magnetization of nickel, by the simple, quantum- 
corrected Weiss theory, with s = 4 and s = 1. Experimental values are those 
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of Weiss and Forrer (15), making use of reductions given by Bozorth (2, p. 720) 
and Myers and Sucksmith (5). The use of slightly different figures quoted by 
Mott and Jones (4) shows that the observed variation of intermolecular field 
is not seriously modified by a different choice of experimental results. Columns 
3 and 4 of Table I give intermolecular field, F, in terms of its value Fy at 6,, 
calculated from the appropriate reductions of the Brillouin function: 


F =| 
aad cee he Sasi“ cs 
[1a] s=}4 tanh( oo? Fe 0; 
adil. 2—y/4—3(e/o0)" _ (.z.. =) _3 
[1b] s=l1: a tanh RT o0/” Fy = 9 ké. 


The first observation from Table I is that the problem is not essentially 
different with s = } and s = 1, although its explanation in terms of ordering 
processes is only possible with s = 1. 

3. NUMBERS OF NEAREST AND NEXT-NEAREST NEIGHBORS 

The face-centered cubic structure may be divided into four completely 
equivalent sublattices. In any unit cell they are represented by the cell corners 
and each opposite pair of cube faces. Each lattice site has 12 nearest neighbors, 
of which there are four on each of the other three sublattices, and six next- 
nearest neighbors, all on its own sublattice. We may find expressions for the 
average number of nearest and next-nearest magnetic neighbors for m magnets 
disposed on four sublattices each with }N sites, and occupied by w, x, y, 2 
magnets respectively, so that 


[2] wtxty+z=n 
The average number of nearest neighbors to each magnet on the w sublattice 
is 


‘WA -+-4 nat wat - . § ty+s). 


Giving this the weight w/m for the w magnets and adding similar expressions 

for magnets on the other sublattices we find the average number of occupied 
£ £ 

neighboring sites for each elementary magnet: 


32 
7 7 (wxt+wyt+wetxy+txst+yz). 
A similar argument gives the average number of next-nearest neighbors: 
, 24 2 2 
¢- 7 (w?+x°+y?+2"). 


These may be combined song equation [2] to give 
n 


ath oR 


We are interested in finding the maximum allowable ranges of variation of 
n and 7’ with values 0.604 and 0.302 for n/N, corresponding to s = } and 
s = 1 respectively in nickel. mmx is realized when w = x = y = 2 = 3m (com- 








664 CANADIAN JOURNAL OF PHYSICS. VOL. 33 


plete disorder) and nmin is found by filling first one sublattice and then another, 
and so on until the required number of elementary magnets have been used 
(ordered state). We first make an assumption, which is reconsidered for the 
calculation in §4, that there is no prece* >‘ ordering within each sublattice, 
but simply a rearrangement of elementary 1 ...' >ts between them. This yields: 





s= j s=]1 

Ordered state: 7 = 7.3 7 =.3.6 
n’ = 3.6 n = 18 

Disordered state: n = 6.1 n = 14 
n = 5.4 n = 5.2 


Although it is a substantial approximation which detracts from the point 
of the present approach, we may consider the intermolecular field at each 
magnet to be that due to interaction with the average number of nearest 
neighbors: 

Fe= nJitn’ J» 


where J; and J2 are nearest and next-nearest neighbor exchange interactions. 

Using the value of F/F», extrapolated to 0°K. and the above values of 7 and 7’ 

we find, assuming the variation of F in Table I to be due to changes in 7 and 7’, 

the following relations between J; and Je: 

s= 4: J2 < —1.04 J; 
I: J2 < —0.2 Ji 


In fact this approximation requires negative next-nearest neighbor inter- 
action, being impossibly large for s = 4. The same difficulty regarding s = 3 
is found in the more exact approach which can be made in the limiting case of 


low temperatures. 
4. THE APPROACH TO ABSOLUTE SATURATION AT 0°K. 


Plotted on a o/oo vs. T/6; curve, experimental results can be fitted to the 
Bloch law (1): 


(equality signs if disorder is complete at 6,). 
sS= 


1—oa/oo = C(T/6)3/? 


up to about 7/@ = 0.3. This fact is remarkable in view of the approximations 
made in its derivation which restrict the range of validity almost to immeasur- 
ably low temperatures, so that the experimental and theoretical laws must be 
regarded as independent. This view has been expressed by Néel (9) as an 
interpretation of the formula developed by Schafroth (11) giving important 
higher powers of (7/@). The same implication can be made from the earlier 
less complete formula of Opechowski (10). Table I gives an alternative method 
of viewing the experimental values which does not indicate a T*/? law, since this 
law would put F = 0 at the 0°K. extrapolation. It thus appears that insen- 
sitivity of the (a/ao) vs. (7/@) curve to the exact function of (7/@) tends to 
mask the nature of their interdependence. 
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The limiting intermolecular field at 0°K. on the four-sublattice model is 
readily calculable and may be compared with that at the Curie point. If 

= 3, then n/N = 0.604 (Bohr magnetons per atom), so that at complete 
order every elementary magnet has at least four others as nearest neighbors as 
compared with the disordered average 7.3. This makes the wide range of 
observed values of F (Table I) impossible to explain without assuming an 
unacceptably large negative value of Jz. Further considerations are therefore 
restricted to coupled spins, s = 1, n/N = 0.302. 

In the ordered state we may regard the w sublattice as being filled, leaving 
0.052N magnets for the x sublattice. Assuming that these are uniformly 
spread throughout the structure, each providing a nearest neighbor for four 
w magnets, none of which has two nearest neighbors, then there are (0.25N 
—4X0.052N) = 0.046N w magnets with no nearest neighbors, experiencing 
only the weak field of next-nearest neighbor interactions. It is reasonable to 
expect that J; >> Jz, whence we may assume that at low temperatures almost 
all of the observed value of (1—a/ao) is accounted for by such magnets, since 
they have the highest probability of being antiparallel to the domain direction. 
Multiplying the values of (1—a/oo) by 0.302/0.046 to find the corresponding 
values for the 0.046N magnets alone and substituting these into equation 
[1b] we can obtain a graphical extrapolation Fy/F,» = 0.08 for these magnets 
alone. Since it is a reasonable approximation to consider the field at the Curie 
temperature to be that due to the average number of nearest neighbors in the 
disordered state we may put 


eee 
3-6): +1-8J_ om, 
whence 
[3] Jo = 0.04 Ji. 


5. SPECIFIC HEAT CONTRIBUTION OF THE SECOND ORDERING PROCESS 


The parameters £;, E, have been introduced without restriction as to mag- 
nitude, which must however be similar to the exchange energies to spread the 
process of ordering of elementary magnets on the lattice sites over the same 
temperature range as the alignment process. The calculations of Néel (8), 
substituting 0, for 6;in 7/0, can be extended to show that this change roughly 
cancels the excess of (F/F») over unity in the range 0.5 < T/@ < 1.0, so that 
to a first approximation we may regard disorder as being complete by about 
T/@ = 0.5. It is however necessary to examine the process to establish the 
absence (at high temperatures) of an ordering ‘‘tail’”’ which could give a specific 
heat contribution and add to (@,—6,). The examination as made is subject to 
Néel’s fluctuation correction, being a consideration of free energy above the 
Curie point. 

Use is made of the requirement that the free energy F of the system must be 
a minimum, where 
[4] F = E-KTInA. 


E is the energy, which if spontaneous magnetization, even locally, is zero, 
reduces to 
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2 
_ , _ 8n',, (8E:—12EF») 

E = 5Q2it1E:) =z A 7 

by neglecting all terms not connected with the ordering. A is the number of 
arrangements of the system with the distribution w, x, y, z of magnets and is 


easily shown to be given by 


InA = -t> ee is in( 42) 


Differentiation of F with respect to changes of the type Ax = —Aw, Ay = 0 
= Az and equating AF to zero for equilibrium at any temperature yields three 
equations connecting w, x, y, z, which with equation [2] may be used to find E. 
Numerical calculations show a \-point discontinuity in dE/dT due to this 
process, and that the system approaches complete disorder catastrophically 
at a critical temperature in a manner similar to the more simple two-sublattice 
model. The peak value of excess specific heat is about 0.8R which specific heat 
measurements (13) show to be far greater than any observed anomaly above 
6, in nickel. It is therefore concluded that disorder is complete at 6,, and that 
E; is of the order 3J,. 

It must be remembered that this approach is not statistical and is subject 
to the Néel fluctuation correction. 


(w°+x?+y?+3") 


6. CONCLUSION 


Tables comparing experimental values of reduced spontaneous magnetiza- 
tion, ¢/ao, as a function of reduced temperature, 7/0, with values calculated 
from the quantum-corrected Weiss theory are misleading in the respect that 
they indicate the important disagreement to be at the higher temperatures. 

Table I gives the values of intermolecular field which are needed to explain 
the observed variation of o/oo in nickel, assuming correctness of the Weiss 
approach in all respects except invariability of the field constant. It shows that 
at 0°K. this constant extrapolates to 1/10th of its Curie point value. This is a 
significant disagreement with the 7*/? law, generally regarded to hold at low 
temperatures, which would require the Weiss field constant to extrapolate to 
zero. The disagreement is not apparent from a direct consideration of the values 
of o/ao. 

The postulated model of wandering elementary magnets, discretely located 
at the lattice sites of a face-centered cubic structure, and subject to a spin- 
independent mutual repulsion as well as a larger exchange interaction, is found 
capable of explaining the 10:1 variation in field constant provided that: 

(i) 0.302 coupled spins per lattice site are assumed, 

(ii) next-nearest neighbors exchange energy is 4% of that for nearest 

neighbors. 

The concept of an ordering process which may occur in pure metals, as a 
consequence of the non-integral Bohr magneton numbers which are observed 
experimentally, is more general than the ordering of alloys as considered by 
Goldman (3) and Smoluchowski. It follows more closely the approach pro- 
posed by Van Vleck and Hurwitz (14), who point out that some measure of 
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compromise with collective electron theory is made necessary by neutron 
diffraction results. 

Such a compromise must come through a modification of the assumption 
that the elementary magnets occupy discrete lattice sites. It simplifies the 
problem to assume instead that they have probability distributions over the 
lattice sites with unequal probabilities of occupying each of the four sub- 
lattices of the face-centered cubic structure. The intermolecular field will then 
be an average over all lattice sites of the field which a discretely located mag- 
net would experience, weighted according to the probability of occupation of 
the sites. This leads to the conclusion of §3 that next-nearest neighbor exchange 
interaction is negative with a value one fifth that for nearest neighbors. 

Calculations for the two sublattices of body-centered iron are simpler, but 
in this case it is difficult to decide how the spins are coupled, being in excess 
of two per atom unless a 10% orbital contribution is allowed, in which case all 
lattice sites are occupied and the simple ordering process visualized here 
cannot appear. 
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LONG RANGE FORCES BETWEEN HYDROGEN MOLECULES! 


By F. R. Britton AND D. T. W. BEAN 


ABSTRACT 


Long range forces between two hydrogen molecules are calculated by using 
methods developed by Massey and Buckingham. Several terms omitted by them 
and a corrected numerical factor greatly change results for the van der Waals 
energy but do not affect their results for the static quadrupole—quadrupole 
energy. By using seven approximate ground state Hz wave functions information 
is obtained regarding the dependence of the van der Waals energy on the choice 
of wave function. The value of this energy averaged over all orientations of the 
molecular axes is found to be approximately —11.0 R~* atomic units, a result in 
close agreement with semiempirical values. 


INTRODUCTION 


Massey and Buckingham (7) have calculated the static quadrupole—quad- 
rupole and the dynamic dipole-dipole (van der Waals) interaction energies 
between two hydrogen molecules, but in making their calculation they omitted 
several terms that have a large effect on the results obtained for the van der 
Waals energy. These terms are included in the present treatment and a numeri- 
cal factor has been corrected. 

The methods used by Massey and Buckingham are outlined in order to 
show where the present treatment differs from the original one and to facilitate 
reading. The original treatment is referred to by numbering in the form 
M.B. 1. 


METHODS 


The Hasse variational method is applied to the function ¥(1+ Av), where 
WV is a wave function for two non-interacting molecules, \ is a parameter, and 
v is the molecular interaction potential. To a second order of approximation 
the minimized interaction energy is given by M.B. 9 as 


2 * _ 2m {fw*o'Wdr— (1/a)( pant) 
[1] _ 1 fw anee h* Wd; (Vv) dr : 


where a = { ¥*W dr and i = 1, 2, 3, 4 are electron coordinates. 
An alternative second order perturbation treatment gives M.B. 10 as 


. * 2 
[2] «= : | f voir 14 J WoWdr—*( Jwtotdr)? + x Elfe.‘eedel . 
where E, is the ground state energy and E, that of another state of the non- 
interacting binary molecular system. 
Electrons 1, 2 are assigned to molecule A and electrons 3, 4 to molecule B. 
The vector joining the mass centers of the molecules is R; d, and dx join the 
nuclei of molecules A and B respectively; r., rz are position vectors of electrons 


1Manuscript received June 22, 19565. 
Contribution from Department of Mathematics, Hamilton College, McMaster University, 
Hamilton, Ontario. 
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in molecule A and f3,r4 position vectors of electrons in molecule B. The 
molecular axes are inclined to R at angles a = (d,, R), 8 = (dg, R) and @ is 
the angle between the projections of d, and dg on a plane perpendicular to R. 


[3] VW = ¥(ri, Fe) ¥(Es, Fs), 


where y is a normalized ground state Hz wave function. d, = dg = do, where 
do is the equilibrium internuclear distance in He. 

Only the terms in the asymptotic expansion of v involving R~® and R“ are 
required, and these are labelled v7; and v2. v; gives the static quadrupole— 
quadrupole energy, varying as R-5, when used in place of v in the first integral 
of [1] or [2]. The van der Waals energy, varying as R~*, is given when v is 
replaced by v2 in those integrals in [1] and [2] that involve v? and in }°,(V; 7)? 
in [1]. The final summation term is neglected in the calculation of van der 
Waals energy from [2]. All other integrals in [1] and [2] are zero to the order 
R-*. 


THE VAN DER WAALS ENERGY* 
lx 
[4] = pe = €p€ 913 (Lp* Ro) (r,"Ro)— (f'f,)], 
a. % 


where p = 1, 2 and g = 3, 4 for electrons; r, = +3d,4,r, = +4ds for nuclei; 
R, = R/R; e, = e, = —e for electrons and e, = e, = e for nuclei. 

Symmetry properties of ground state Hz wave functions are used after 
squaring [4] and after squaring V; v2 to give expressions for use in [1] and [2] 
as follows: 

: a : 
5] Ge = “pil tats)? (Rata) (tata) + (Pres) (Fes) + (rts) (Fata) 
+9{ (ri a Ro)*(r3° Ro)°+ (r, ‘ Ro) (fo: Ro) (rs: R,)° 
+ (£3: Ro) (t4-Ro) (ti Ro)’+ (r1-Ro) (t2° Ro) (fs Ro) (f4" Ro) } 
—6{ (fi-fs) (f1- Ro) (ts-Ro)+ (f1-8s) (f1- Ro) (ta Ro) 
+ (fi-T3) (Te: Ro) (r3-Ro)+ (fi-fs) (fe- Ro) (rs-Ro)}], 


and 


4 
[6] (7.02)? = FEol2rn?+2 (e142) +3 (G1-Ro)*+3(ts-Ro)* +3 (Fx: Re) (Fe Ro) 


i 


+3 (r3° Ro) (rq: R,)]. 


If molecular symmetry properties are used before squaring [4] and V; 2, 


*The form [4] is used for ve in order to compare our revised results with those given by Massey and 
Buckingham. The development is made somewhat easier by using the Cartesian form 


2 
"5 x & Cplq(2prq—XpXq—VpVe)s 
q 


where the z axes are in the direction R, and transforming coordinates for vz by rotation of the 2 axes 
to axes along the molecular axes. 








670 CANADIAN JOURNAL OF PHYSICS. VOL. 33 


only those terms given by M.B. 14 remain, three terms of [5] and three terms 
instead of the six in [6]. 
Using the notation 


(fa = Sftvan, r2)]'f dridre, 
(fs = SS (vrs, r.)]'f dr3dr4, 
(f)) = SSSP lon 2) Wes, ra) Pf dridradradrs, 


and cosine formulae of spherical trigonometry it can be shown that 
(ri-T2), = (f1-T2 cos(fi, dy) cos(f2,d,4)),4 = X, 
(ri?)a = (ri? cos?(f1, da)) a+ (ri? sin?(f1, da)), = C+S, 
((t1-Ro) (f2-Ro)), = X cos? a, 
((13-Ro) (f4-Ro))s = X cos? B, 
((f1-Ro)*), = $S+(C— }$S)cos? a, 
((ts-Ro)?)p = 35+ (C— }S)cos? B, 
[7] (((t1-13)?)) = (C— 4S)? cos?(da, dg) +S(C+4S), 
(((f1-83) (F2-F4))) = X* cos*(dy, ds), 
(((fi-Fs) (2-s))) = X?{ (C—4S)cos?(d,,dp)+ 35}, 
(((f1-Ts) (Pi Ro) (f3-Ro))) = 4S?+ 35(C— 3S) (cos? a+cos? B) 
+(C— 4S)? cos a cos 8 cos(d,4, ds), 
(((f1-%3) (Pi: Ro) (f4-Ro))) = X {4S cos? B+ (C— }$S)cos a cos 8 cos(d,4,dz)}, 
{{(f1-13) (2 Ro) (ts-Ro))) = X {45 cos? a+ (C— }$S)cos a cos 8 cos(da, ds) }, 
(((f1-83) (fe: Ro) (f4-Ro))) = X* cosa cos B cos(d,, ds). 


From [5], [6], and [7] 


[8] (((v2)*)) = “alS(X+C+S) +4 S(X+C—45) (cos’a+cos’) 


+(X+C—45)*(sin a sin B cos @ — 2 cosa cos B)’], 
and 


9] (Xe (W.02)*)) = Gl(X+C-45){2+3(cos'a + cos"8)} +65]. 


The integral X does not appear in the Massey and Buckingham treatment, 
and the factor 3 in [9] replaces a factor 6 in M.B. 23. 

From [1], [2], [8], and [9] the van der Waals and quadrupole—quadrupole 
energies are given by 


_ _2( ao )* eal fla, 8, 0)" fay, 
= —£(#) h(a, B) To, R 2 8(a, B, 8), 


[10] 
e(a\*, e'( a )° ’ 
€p = i y Ci fla, B, jat(es) C2 g(a, B, 6), 


where 
f(a, B, 0) = 14+ M (cos? a+cos? 8)+ N(sin a sin B cos @—2cos a cos B)?, 
h(a, B) = 1+ L (cos? a+cos? B), 
g(a, 8,0) = 1+17cos* a cos? B—4sin 2a sin 28 cos 6—5 (cos? a+cos? B) 
+2sin? a sin? B cos? 6, 
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_ 3(X+C—3S) ue = 3X4 0-3S) ae (X+C-45)* 
2X+2C+5S ’ 2(X+C4+S) ’ S(X+C+4+S) ’ 
_ 8S°(X+C+5)’ , _ 4S(X+C+5S) ala adele 
c= 2X+2C+5S ’ Gy = Eo ’ —_ 3(C 35 ad ) e 


Compared with M.B. 33, which corresponds to [10], the functions g(a, 8, 6) 
and cz are unaltered. They define the R~ static quadrupole—quadrupole energy 
and details of their derivation are given by Massey and Buckingham. The 
large sinusoidal variation of this energy with orientation of the molecular axes 
is illustrated by M.B. Figs. 2A and 2B. Averaged over all orientations this 
energy vanishes so that for gas kinetic phenomena in which symmetry prop- 
erties permit this energy to be neglected the van der Waals R-* dynamic 
dipole-dipole energy is of predominant importance provided that R is suffi- 
ciently large to neglect the R-* dynamic dipole-quadrupole energy. 





Fics. 1 and 2. Variation of f(a, 8, @) with orientation of molecular axes as given by the 
Coulson (C) and Wang (W) functions compared with that given by M.B. Figs 2A and 2B 
which is shown by broken curves. (Note that the labelling, @ = 0° and 6 = 90°, should be 
interchanged in M.B. Fig. 2B.) 

Massey and Buckingham used the Coulson (1) and Wang (9) functions for 
Hz in their calculations. In [10] as compared with M.B. 33 the function X¥+C 
replaces C in all constants except ¢2, and L is reduced by a factor 4. The effect 
of these revisions on values of f(a, 8, 8) is illustrated by Figs. 1 and 2 in which 
Massey and Buckingham values are shown by dotted curves. X = 0 for the 
Coulson function and the differences in values of f(a, 8, 6) are due to our cal- 
culations being made for values of constants given by Coulson instead of for 
- those used by Massey and Buckingham which are not carried to as many figures. 
X 0 for the Wang function and our values of f(a, 8, 8) are much lower and 
vary much less with orientation of molecular axes. 

For the Coulson function the additional factor 4 in L does not change the 
character of the variation of the van der Waals energy, Ey4.y,., with orientation 
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but its value averaged over all orientations, Eyqw., as given by e, in [10] is 
— 15.7 R-‘ instead of — 15.0 R-* as given by Massey and Buckingham. (Atomic 
units are used here and following.) 

For the Wang function our revisions produce large changes. Ey.aw, varies 
very little with orientation and e¢, gives Eyaw. = —11.3 R-® instead of 
—15.9 R-°, the additional factor of } in L that tends to increase the magnitude 
of Eyaw. being more than compensated for by X ¥ 0. 

Because the Coulson and Wang functions give greatly differing values of 
E,.a.w. and differing amounts of variation of E,.4.~. with orientation, calcula- 
tions were made for five other functions in order to investigate the dependence 
of van der Waals energy on the function used to approximate the ground state 
energy of He (i.e., }£o). The Heitler-London (3), Gurnee—Magee (2), Rosen 
(8), and the two Weinbaum (10) functions were used. All seven functions are 
given in an appendix together with formulae for the constants X, C, and S 
that appear in [10]. Numerical values of these constants and of ¢; and c,;’ are 
given in Table I. Atomic units are used and functions are referred to in tables 
and graphs as follows: Heitler-London (HL), Wang (W), Gurnee—Magee (GM), 
Rosen (R), Weinbaum Wang-ionic (WWI), Weinbaum Rosen-ionic (WRI), 
and Coulson (C). 



















TABLE I 
Functions HL W GM R WWI WRI Cc 
a —0.3793 —0.3360 —0.2764 -—0.1768 —0.2097 —0.1492 ©) 
(0) 0 
S 2.0804 1.5341 1.4971 1.5386 1.4789 1.4706 1.5744 
(1.534) (1.594) 
X+C-—3S 0.0393 0.0316 0.0286 0.0235 0.1439 0.1050 0.2534 
(0.364) (0.200) 
L 0.0094 0.0102 0.0095 0.0076 0.0471 0.0349 0.0764 
(0.220) (0.120) 
M 0.0187 0.0203 0.0187 0.0151 0.0914 0.0682 0.145 
(0.203) (0.116) 
N 0.00024 0.00035 0.00024 0.00015 0.0059 0.0032 0.0156 
(0.032) (0.0097) 
C1 27.52 11.055 10.355 11.095 10.656 10.228 13.62 
(13.7) (13.7) 
6.25 6.09 5.90 7.26 











; : . 5.96 





M. B. values are given in parentheses. 





Values of Ey.4.w. given by the variational method are more reliable than those 
given by the perturbation method because of the neglect of summation terms 
in [2], and variational values will be considered hereafter. 

The dependence of Fy.4.w. on the relative orientation of the molecular axes 
is illustrated by Fig. 3 and again by Figs. 4 and 5 from which the Heitler- 
London and Coulson results are excluded. Figs. 3 and 4 are for the orientation 
that gives extreme variation of Ey.q.y.,a@ = 8,8 = 0. The perturbation method 
gives similar but smaller dependence on orientation. The values of Ey.a.w. 
given in Table II were calculated analytically and range from —28.1 R-® for 
the Heitler-London function to —10.5 R-* for the Gurnee—Magee function. 
A good approximate value of E,.4.w. is obtained by evaluating c, f(47, 41, 0). 
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Fic. 3. Variation of Ey.q.w. with orientation of the molecular axes for a = 8, 6 = 0. Units 
are atomic. 


= ww 
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Fics. 4 and 5. Variation of Ey.a.w~. with orientation of the molecular axes for a = 8, 0 = 0 
and for 8 = a+47z, 9 = 0. Units are atomic. 


DISCUSSION OF RESULTS 


Hartree atomic units are used in the following discussion. 

The Heitler-London, Wang, and Gurnee—Magee functions for the ground 
state of a He molecule are of the same form (see Appendix). In the case of the 
first two of these functions the 1s orbitals are centered at the two molecular 
nuclei but in the case of the latter function the 1s orbitals are centered at two 
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points on the molecular axis each of which is offset a distance x from a nucleus 
towards the center of the molecule. Data available for the Gurnee—-Magee 
function did not permit exact determination of x, internuclear separation d, 
or screening parameter z for the position of minimum energy. Calculations 
were made with small variations from the approximate values x» = 0.065, 
dy = 1.425, 2) = 1.179 and gave small changes in E,.4.y. with values numerically 
smaller than those given by the Wang function. 

Each of these functions gives very little variation of E,.4.~. with orientation 
(see Fig. 3) but the Heitler-London function gives a much larger value of 
|\Zy.aw.| than the other functions (see Table II). However, Hirschfelder and 
Kincaid (4) have shown that the Heitler-London function does not satisfy the 
virial theorem of quantum mechanics and so values of Ey.4.~, given by it cannot 
be reliable. 

















TABLE II 
Ey.a w 
Function do Zo Eo _—_— —_—_—_—_—_—— 
Pert. Var. 
HL 1.518 : 2.232 —11.78 —28.1 
W 1.406 1.166 2.278 —6.25 —11.3 
— (8.11) — (15.9) 
GM 1.425 1.179 2.308 —6.04 —10.5 
R 1.416 1.170 2.297 —6.31 —11.2 
WwIl 1.417 1.193 2.296 —6.48 —11.7 
WRI 1.416 1.190 2.303 —6.18 —11.0 
C 1.40 2.266 —8.03 —15.7 
— (7.62) — (15.0) 
Observed 1.40 2.351 


The Wang function can be considered to be of Gurnee—Magee type with 
x = 0, and both functions satisfy the virial theorem. For such functions 
Ey.a.w. depends on Zo and the equilibrium position separation of orbital centers 
(see the appendix and [10]), and, by using data given by Gurnee and Magee, 
it can be shown that |Fy,4.w.| has a maximum value for an orbital center separa- 
tion of about 1.40 at the equilibrium position. For the Wang function the 
equilibrium center separation is also do and is slightly larger than 1.40, so that 
this function gives a value of |Eyay.| very near to the maximum given by 
functions of Gurnee-Magee—Wang type that satisfy the virial theorem. This 
value is 11.3 R-®. 

The Wang function gives a value of dy only a little larger than the observed 
value but takes no account of polarization and ionic effects which tend to 
increase the concentration of electron density towards the center of the 
molecule and therefore to reduce values of |Ey.aw.|. The effect of adding terms 
that approximate polarization is seen in results obtained for the Rosen func- 
tion. The greater concentration of electron density towards the center of the 
molecule produces less variation of Ey.4,w, With orientation (see Fig. 4). Both 
do and 2o are larger than for the Wang function and whilst the increase in do 
tends to increase values of Ey, the increase in zo has the opposite effect. 
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Ev.aw.| = 11.2 R-S, a little smaller than the value given by the Wang func- 
tion. The effect of adding weighted ionic terms, as in the Weinbaum Wang- 
ionic function, is to produce greatly increased variation of Ey... with orienta- 
tion (see Fig. 4). Both zo and dp are larger than for the Rosen or Wang functions 
and |Eyaw.| = 11.7 R-*&. The Weinbaum Rosen-ionic function includes both 
polarization and less heavily weighted ionic terms. dp is the same as for the 
Rosen function but 2 is a little larger. The variation of Ey... with orientation 
is smaller than that given by the Wang-ionic function (see Fig. 4) and 
|Evaw.| = 11.0 R-*. 

The molecular orbital Coulson function gives greater variation of Eyay. 
with orientation than any other function considered (see Fig. 3). |Evaw.| 
= 15.7 R-®. This function uses a value of dy that differs very little from the 
observed value and the constant used in exponential factors (see Appendix) 
is derived from the James and Coolidge function (5). However, the Coulson 
function does not take account of the electron interaction terms that were 
shown to be so important by James and Coolidge. Molecular orbital functions 
of the form ¢(1) ¢(2) give X = 0. They are not well suited to represent ground 
state functions of homonuclear molecules and the poor value of Ey obtained by 
using an almost exact value of dy and four constants determined by the self- 
consistent field method indicates a poor representation of electron density. 
Values of Ey.ay. given by this function cannot be reliable. 

The results given by the Coulson, Heitler-—London, and Weinbaum Wang- 
ionic functions will not be considered further, the latter results being dismissed 
because the Rosen-ionic function can be expected to give a better approxima- 
tion to the true wave function. 

Of the remaining functions, that of Wang gives the best dp and poorest Eo 
whilst that of Gurnee—Magee gives the poorest dy and best Eo. The former has 
insufficient electron density towards the center of the molecule and the latter 
probably too much. Therefore |Ey.aw.| probably lies between 11.3 R-® and 
10:5 R-“. 

The Weinbaum Rosen-ionic function gives Eo slightly smaller than the 
Gurnee—Magee and a much improved dy = 1.416. This function can be ex- 
pected to give the best approximation to the true electron density distribution 
of any of the simple functions considered and hence the best description of the 
variation of Ey.aw. With orientation and the best value of Ey... The maximum 
variation from the mean is about 10% which is much less than that given by 
the Coulson function and also much less than that given by Massey and 
Buckingham. The one-parameter Wang function appears to give a fairly good 
value of Eya., but insufficient variation of Ey... with orientation. 

From the above considerations it appears that Eyaw. = —11.0 R-® approxi- 
mately. 





COMPARISON WITH EXPERIMENT 


The most reliable semiempirical calculation of Ey... is probably that made 
by Margenau (6) who based his calculation on oscillator strengths and their 
associated energies that had been determined empirically by Wolf and Herz- 
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feld (11), who fitted the experimental index of refraction of Hz as a function 
of frequency very accurately by a two-term dispersion formula. Margenau’s 
two-term value is Ey4y. = —10.9 R-*, which is in good agreement with our 
value for Eyaw.- 
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APPENDIX 
Hz. WAVE FUNCTIONS AND FORMULAE FOR CONSTANTS 

Molecular nuclei are labelled A, B and electrons 1, 2. 7,4;, 73; are vectors to 
4 = 1 or 2 from A and B respectively. The equilibrium separation of nuclei is 
d. Prolate spheroidal coordinates are defined by 

hi = (raitrsi)/d, wi = (rai—rsd/d, o; = azimuthal angle, 


and were used to derive formulae for constants. Atomic units are used with 
normalized 1s hydrogen-like orbitals of the form 
* es 1 
a(t) = (23/mr)? exp(—2rai). 


The following functions appear in constants: 
oe 
A = (itstis*)e-, 
/ ee ee ee eee 
P = (1+s+3s +13 Je", 
9 2 7 3 1 4) —s 
Q (1+s+26s° +605 +e05 Je’, 
: 9) 08. 28 8 8 8 Be 
G = (1+s+i65 "+1268 +2105 — 9455 — T3905 Je’, 
45.2, 17.3, 19.4, 19 5, 1 6) — 
H = (l+s+us +s +3365 +1605 +7205 )e*, 
9 
Q—P = As’ /20,s = 2d. 


The Heitler-London, Wang, and Gurnee—Magee Functions (8, 9, and 2) 
y(1, 2) = N{a(1)b(2)+a(2)d(1) }, 
N-* = 2(1+-A?), 
X = —N’s?/22?, 
C = (2N?/z?)(14+4s?+A0), 
S = (4N?/z*)(1+AP), 
X+C—43S = N?A?*s?/102?. 


For the Heitler-—London function z = 1, d = 1.518. 
For the Wang function z = 1.166, d = 1.406. 








BRITTON AND BEAN: LONG RANGE FORCES 677 


For the Gurnee—Magee function z = 1.179, d = 1.295 were used, with d the 
distance between the centers of 1s orbitals for an equilibrium internuclear 
distance 1.425 and x = 0.065 in the case of this function. 


The Rosen Function (8) 
¥(1, 2) = N{U,(1) Up(2)+ U,(2)Uz(1)}, 


where 
U,(t) = a(t) {1+ 3dw(1+A; u:)}, 
Us (1) = b(t) {1+ 3dw(1—d,; ws) }, 
s= 1.170 d= 1.4146. c= lz. 


2\ 2 272 2 
Nn 24 (148%) +FY, Fe (1.004% )o2 P, 


d #) 20? 
= ONS a Pas re 
2N ‘(14% "ae 


2 


2 272 2 2 
(imei), fasano-Se}] 


._ 4N? yi a) {( =f | 1 22, 30° \] 
S=-3 [ (1455 1+52 +F\ I+> Ptgwd Q-a2H : 


ll 


>< 
| 


ll 
| bo 
a> 
iw 
, 
_ 
+ 
tJ ie 





2 
5S 2N’ of deo 
#+€-] = a seater” 3 
wd\*As* 3 32, 3 \] 
+7 (1484) "48 g 2d Q—52(8G—H) ; 


The Weinbaum Functions (10) 
(i) (1,2) = Nwtk yi), k = 1/3.9, d = 1.417, where pw is the un- 
normalized Wang function and 

Wr = a(1)a(2)+0(1)b(2). 

z = 1.193 for all 1s orbitals. 
N~? = 2{(1+k?)(1+A?)+4kA}, 

Ps = N?{Xwt22kXwitk?X 7}, 
where Xw, Cw, and Sw are the X, C, and S functions as given for the Wang 
function without the normalizing factor N?. C,S, and X+C—4S are given 


similarly. 

Xwi — 0, 

Cwr = (2/2?) {(1+4s?)A+Q}, 

Sw1 = (4/2")(A+P), 

(X+C—4S)w1 = 3As?/52?, 

X, = s3/23?, 
Cy = (2/2?) (1+4s?+ AQ), 
S; = (4/2)(1+AP), 


(X+C—3S)1 = (s?/2) (1+ 754"). 
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(ii) Y(1,2) = N@eatk yi), k = 1/5.7, d = 1.416, » = 0.0833, where yp is 
the unnormalized Rosen function and z = 1.190 for all 1s orbitals. 


2\ 2 2 , 
rod (os) f(s) sce svansyeaaod)]} 


X, C, and S are as given for (i) with the Wang function replaced by the Rosen 
function and 


Xxx = (2w/2*) (d—2w/2”)Q, 
Car = (2/2?) (1+ gwd) { (1+ 4s°—wd)A+Q}, 
Sri = (4/2*)(1+ 30d) (4+P), 
(X+C—4S)r1 (2/2?) { (1+ gwd) (;55°—wd) A+ (wd — 2w?/2?)Q}. 


The Coulson Function (1) 


where 
$(1) = 0.87758e—°7™: (1+4+-0.27787 1 — 0.12863A,+0.012503A 1? — 0.039589) 11”). 


The following reduction formulae were used in the calculation of constants, 
numerical values of which are aie in Table I: 


— 


SOON pO 


1 
Jiu du 8 n even 


J oh = i ote EON TAN. 
1 
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FIELD THEORETICAL DESCRIPTION OF THE 
SUPERCONDUCTING STATE! 


By Rosert E. PuGH 


ABSTRACT 


The superconducting state is described in terms of the electronic plasma in 
interaction with the electromagnetic potentials. The interaction of the plasma 
with itself and with the lattice is represented by a potential energy determined 
such that the equations of motion are self-consistent. The resulting non-linear 
field equations are solved in the linear approximation. There exist four indepen- 
dent modes of vibration of the field variables corresponding to the plasma oscilla- 
tions, and one longitudinal and two transverse oscillations of the electromagnetic 
potentials. The contribution of these modes to the specific heat is discussed. 


I. INTRODUCTION 


In this work an attempt is made to describe the phenomenon of supercon- 
ductivity by representing the so-called ‘‘superconducting electrons’ as a 
classical charged fluid. It will be shown that the hydrodynamic equations 
describing the fluid are consistent with London’s phenomenological equations 
(3). Since as yet the behavior of the specific heat of a superconducting metal 
below its transition temperature has not been satisfactorily explained, we shall 
try to describe this phenomenon in terms of the energies of the wave motions 
of the electric and magnetic field potentials. We shall begin our discussion by 
reviewing briefly some of the experimental facts associated with the phenome- 
non of superconductivity. 


A. The Experimental Facts 

In 1911, Kamerlingh Onnes (5) discovered that at a certain very low tem- 
perature, 7,, the resistance of certain conductors dropped sharply to zero, 
indicating that no electric field could be maintained inside this ‘superconduc- 
tor’. Kamerlingh Onnes further discovered that this phenomenon of supercon- 
ductivity was destroyed by the application of a sufficiently strong magnetic 
field, H,(T). 

Considering the superconductor to be a perfect conductor (infinite conduc- 
tivity), then it follows from Maxwell’s equations that the magnetic field inside 
a superconductor is constant (B = 0). One would then expect that if a super- 
conducting metal were placed in a magnetic field and then cooled below its 
transition temperature, the magnetic field would be ‘frozen in’. It was found, 
however, by Meissner (4) in 1933 that when the metal became superconducting 
the field in the neighborhood of the specimen increased, thus showing that the 
field inside the superconductor was pushed out, leaving no magnetic field 
inside. It was also shown that when a superconductor initially in zero field is 
placed in a magnetic field, the field inside the superconductor remains zero. 
Hence one concludes that this process is a reversible one and that the magnetic 
field may therefore be treated as a thermodynamic variable. 

1 Manuscript received May 19, 1955. 


Contribution from the Department of Physics, University of British Columbia, Vancouver, 
British Columbia. The author is a holder of a National Research Council Bursary, 1954-55. 
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Treating the magnetic field H as a thermodynamic variable of the state of 
the superconductor, one can easily derive from the first and second laws of 
thermodynamics that the difference between the specific heat in the super- 
conducting state and the specific heat in the normal state is 


C.— eC a 
Tv 


The specific heat at constant volume in the normal state at low temperature 


can be represented by 
C, = aT+464.4(T/6)'. 


The T* term, which is hereafter denoted by Cz, is attributed to the lattice 
vibrations according to the well-known Debye theory of specific heats. The 
linear term has been shown (6) to result from the application of Fermi—Dirac 
statistics to a free electron gas, and is hence usually referred to as the normal 
electronic specific heat and is designated by C,". 

It is generally found that at the transition temperature 7, where the change 
from normal to superconducting states takes place, the specific heat as a 
function of temperature is discontinuous (see Fig. 1). It is seen from Fig. 1 that 


Q__ MOLAR HEAT 





01 


Fic. 1. Molar heat capacities of niobium. 


the specific heat in the superconducting state, C,, is greater than the lattice 
contribution, C,, at all points, but not always greater than the normal elec- 
tronic specific heat. Hence it is generally assumed that the lattice vibrations 
still contribute to the specific heat in the superconducting state, whereas the 
normally conducting electrons do not contribute. Since the difference C,— Cz, 
cannot be accounted for by C,°!, this difference is called the superconducting 
electronic specific heat, and is designated by C,*. Experimental measurements 
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of C,*' generaily indicate that it follows a T* law. Recent work on the specific 
heat of superconducting niobium (1) shows that in this case C,°' rises much 
faster than any 7* dependence. We shall later discuss the contribution of 
various modes of vibration to the superconducting specific heat, and will give 
a possible reason for this deviation from a 7* dependence. 


B. The Present State of the Theory of Superconductivity 

There exists at present no atomistic theory of superconductivity. One can, 
however, specify some of the conditions that a successful theory of supercon- 
ductivity must fulfill. It must 

(i) agree with the phenomenological theory (London’s equations), 
(ii) give the transition temperature and the critical magnetic field, and 

(iii) give the isotope effect (M? T. = const.). 

Taking the usual view that superconductivity is caused by the conduction 
electrons, and is not a property of the lattice, one can say that a successful 
approach towards the problem has been made if one can explain the origin of a 
large gap in the energy levels of the conduction electrons near the surface of 
the Fermi sphere. In fact, if electrons were Bose particles, one could construct 
a fairly reasonable model of the superconductor, in which the existence of the 
critical magnetic field is explained in terms of the Zeeman effect field which 
would tend to close the gap between the ground state and the first excited state 
of the electrons in a superconductor of reasonable dimensions. This effect 
comes out in the right order of magnitude in the approximation in which the 
interaction between the electrons themselves and between the electrons and 
the lattice is neglected. Such a model is unfortunately of little use in the case 
of a Fermi gas of electrons, since the level density near the surface of the Fermi 
sphere is so large that the magnetic field needed to make neighboring energy 
levels coincide by Zeeman effect is extremely small (10~* oersteds). It can 
therefore be only by taking into account the interaction between electrons 
themselves, and between the electrons and lattice vibrations, that one can hope 
to obtain a successful start towards an atomistic theory of superconductivity. 
The very existence of the isotope effect indicates that the width of the energy 
gap at the surface of the Fermi sphere depends essentially on the interactions 
between the electrons and the lattice vibrations. 

It is beyond the scope of the present work to try to describe the mechanism 
which sets up the superconducting state, which is presumably due to the inter- 
action of electrons among themselves and with the lattice ions. It will be 
assumed that this interaction leads to a rigid ground state for the electrons, 
i.e., most of the electrons must be in their lowest states, and the probability 
that they can rise to states of higher energy must be small. 

In the following, the electronic description of the superconductor will be 
characterized by a system of continuous electric charges which interact with 
the electromagnetic field which they themselves produce. In other words, the 
electronic behavior is characterized by a system of simultaneous differential 
equations consisting of Maxwell’s equations and the classical hydrodynamic 
equations for a continuous fluid. 
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II. THE LAGRANGIAN FORMULATION OF THE SUPERCONDUCTING STATE 
AND THE RESULTING EQUATIONS OF MOTION 


We shall try the following Lagrangian density :* 
° 1 Lie: > * 
o 2 a A 2 2 © . 
[1] L = —po—ipv ter (cE —H ee $+V ) —— pi 


where H is the magnetic field, E the electric field, @ and A the electric and 
magnetic vector potentials respectively, « the dielectric constant, and Vv 
(defined by equation [2] below) the velocity of the fluid. ¢ is the velocity 
potential, and p is the mass density of the fluid. The first two terms are the 
usual classical hydrodynamic terms for a free continuous fluid of density p; 
by themselves, they give upon independent variation with respect to p and ¢ 
the equation of continuity and the Bernoulli equation for a continuous fluid. 
The last three terms are the electromagnetic terms which give, by themselves, 
upon variation with respect to A and ®, Maxwell’s equations. There is, how- 
ever, one important difference: the usual term (e/m)p® has been replaced by 
(e/m)pi%;, where p; and 4, are the deviations from the average density po 
and the average potential &) respectively. This means that the background of 
positive lattice charge effectively cancels out the Coulomb interaction of the 
fluid leaving only the interaction between the net electric charge density 
(e/m)p:, and the net electric field potential 4). p and ® are imagined to be 
written as p = potpiand ® = &)+ 4). 

In addition we define the following: 


; ss cee 
[2] v= Vo mc A, 
[3] H = curl A, 
[4] E = -va—*A. 


The choice [2] has been made to satisfy identically London’s equation 


[2a] curl Ajz = 1H where Ajz = =e, 


while [3] and [4] are the usual expressions for the electric and magnetic fields 
in terms of the potentials ® and A. 

We will now assume that the density of the fluid is almost constant and that 
variations in density are small, i.e., we assume that in the expansion p = po+p1, 
the deviation p; is much smaller than po, and that po is constant. This assump- 
tion, with the previous assumption that the positive lattice charge effectively 
cancels out the electric field due to the electrons, makes it reasonable to assume 
that ¢, A, and 4; are small, so that second order terms in them can be 
neglected. 

*This Lagrangian is chosen for simplicity; one can try another Lagrangian which is derived 
from analogy with the classical Schroedinger field, but tt turns out that this more complicated 
type of Lagrangian leads to the same solutions if one demands that London's first equation be 


satisfied. The treatment of this type of Lagrangian is outlined in Appendix A. 
{This is shown to be true in Appendix A. 








: 
‘ 
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By independent variation of the Lagrangian density with respect to ¢, p, ®, 
and A, we get the following equations of motion: 


[5] pit eov"e——— poV-A = 0, 

(6] g+< %, = 0, 

” 2, € 4 , 47 ay 

[7] Vv A-—aAt mc n(ve mc A) = 4%, 
2», —£ 4 tre 

[8] V'b— 4 b+ p, = 0, 


neglecting second order terms. Equation [5] is the usual equation of continuity 
expressing the law of conservation of matter. Equation [6] is Bernoulli's 
equation for the conservation of energy in the case that the kinetic energy of 
the fluid is negligible. Equations [7] and [8] are the usual electromagnetic 
field equations. London’s first equation follows directly from equation [6]. 
We have 


| ms 


[6a] dAjez _ 9(m/e)v _ my = 


: 1: 
at at A= “88-74 =E 


5 


which is London’s first equation. 
The equations [5] to [8] may be solved for p; as follows: differentiate [5] 
with respect to time and substitute ¢ from [6]. Making use of the Lorentz 


condition V-A = —(e€/c), equation [5] becomes: 

4ne° 
9 fi: +~— pop: = 0. 
[9] pit me POP2 0 
Hence , pi = f(r) exp (twot) 
where 

2 4ne* 

[10] wo = oF Po 


and f(r) is, as far as equation [9] is concerned, an arbitrary function of the 
co-ordinates. However, in order to solve the rest of the equations, we must 
know the explicit form of f(r). It turns out that if f(r) is taken to be a constant, 
the Lorentz condition cannot be satisfied. In Appendix A, where a more general 
Lagrangian is considered, f(r) is shown to be of the form exp(iKy-r) where 
Ky is a vector whose magnitude depends on the coefficient of a term which 
does not appear in the Lagrangian used above. It is therefore convenient to 
take 

[11] f(r) = const. exp(iKo-r) 


where for this discussion Ky can be taken to be any arbitrary non-zero vector. 
Hence we have: 
[12] pi = const. exp(1Ko-r+iwof). 
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We may now solve equation [8] for #,. We get: 


, ; 4 re 
[13] , = const. exp(iK - Ebel) +E S— elo en] pi 
where 
[14] K? = ew?/c?. 


For brevity we put 9%,’ = const. exp(iK-r+iwt) and a = 4me/[me{Ko?— 
(€/c?)wo?}]. Equations [5] and [6] yield the following solution for ¢: 

e ee ee ee , 
[15] o= cont tt weal o;+¢ 
where ¢’ is a general solution of ¢’ = 0 and V°¢’ = 0. Finally equation [7] 
yields: 





[16] A = Vd +46: +s Véi-+Ao exp(iKy r+iwst) 
where Ap is a constant vector, and 

eK? 
[17] 3” = . 1 en”. 


It is seen from [13] and [16] that the Lorentz condition is satisfied provided 
that Ao-K, = 0, i.e. provided that the vector potential 
[18] A’ = Ay exp (1K, -r-+ iw!) 


is a transverse wave. 
It should be pointed out that the magnetic vector potential A given by 
[16] has a very unusual feature: it has a longitudinal component, (c/w*)V%y1'. 
We see from the above that the solutions are constructed from four indepen- 
dent modes of oscillation: p:,-#,', and the two components of A’ which are 
orthogonal to its propagation vector K;. p; has a fixed frequency, wo, given by 
[10], and an arbitrary propagation vector, Ko. 4,’ has a frequency, w, and a 
propagation vector, K, which are related by 
w? = c?K?/e. 


A’ has a frequency w, and a propagation vector Ky, related by 
wi? = (c?K,?/€)+ wo’. 


It is shown in Section III that, for thermodynamic considerations, wo? is 
negligible compared with c?K,?/e. The contribution of these four modes of 
oscillation to the specific heat of the system is considered below. 

Equation [10] is the clue to a possible explanation of the high specific heats 
of superconductors. It states simply that the frequency of the classical charged 
fluid representing the superconducting electrons, i.e., the frequency of the 
so-called plasma oscillations, is a constant. We could then imagine the 
following model. In the normal state, the conducting electrons are essentially 
free, as shown by the existence of C,°'. At some critical temperature the elec- 
trons, under the influence of the lattice ions, vibrate with some fixed frequency 








PUGH: SUPERCONDUCTING STATE 685 


which is a subharmonic of the maximum lattice frequency, i.e. a subharmonic 
of the upper limit of the Debye spectrum of lattice frequencies. Let us assume 
that this plasma frequency is of order of magnitude kT ./h ~ 10! to 10". 
This frequency results in an extremely large dielectric constant as was pointed 
out by Ginsburg (2). For example, if we take the ‘‘number of superconducting 
electrons” per unit volume po/m to be 5X10?!, then from equation [10] we 
have 


4ne* 8 10 
€=-z773po~ 10 told. 
Mm Wo 


This large dielectric constant plays a very important role in the behavior of 
the superconductor. Such a large ¢€ means a very small velocity of light 
(c/ Ve ~ 10° cm./sec.). Now the smaller the velocity of light, the larger is 
the momentum associated with black body radiation of given energy E. 
Hence, since the number of photons which can have energy between E and 
E+4dE is proportional to p*dp, there are more photons with energy E when e 
is large than when ¢ is small. In fact, the number of photons of given energy is 
proportional to e*/?. Because of the greatly increased number of photons 
present, a much greater expenditure of energy will be required to raise the 
temperature of the superconductor. This means a larger specific heat in the 
superconducting state than in the normal state. This subject is treated 
quantitatively in the next section. 


III. SPECIFIC HEAT 


From the solutions for the equations of motion given in Section II, we see 
that there are four independent modes of vibration in the superconductor. 
They are: 

(1) the plasma oscillation, 

(2) the electric field potential oscillation #,’, 

(3) and (4) the two independent modes of vibration due to A’. There are 
two modes because A’ is a transverse wave and hence there are two 
independent mutually orthogonal oscillations each perpendicular to the 
direction of propagation of the wave. 

In thermal equilibrium, each of these modes will contribute separately to the 
total energy of the system. We will first show that the contribution of the 
plasma oscillations to the specific heat of the system is negligible and then 
proceed to calculate the contribution of the other three modes of vibration. 


(1) Plasma Contribution 
It was shown in Section II that the density of the continuous charged fluid 
representing the superconducting electrons was of the form 
P = potpi = potconst. exp (7Ko-r+-iwol). 
We define the energy of such an oscillation as Ey) = hwy and the momentum as 
P, = hKo. Then the total energy per unit volume due to the plasma oscillation 


is 


: 4x?" EyPo dPo 
Foor = Qah) Jo geht 4 


(7 
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where Pomax is the maximum momentum associated with the plasma oscilla- 
tions, k is Boltzmann’s constant, and T is the temperature. Because Eo is 
independent of Po, the size of Eoto, will be determined entirely by Ey and 
Pomax- Since we know nothing about the size of Pomax, the best we can do is to 
assume some upper limit for it. We will assume that 


“ 4 
Pomax = ws Eo = a ten) ; 


c c\ m 
presumably it would be much smaller than this. We therefore have: 


Eotot = AP mask o : 
3(2ah) (e°’” —1) 
We will now assume for the sake of the following order of magnitude calcula- 
tion that Pomax and EK» are temperature independent. In general, of course, 
they will not be independent of temperature, but taking this into account would 
not change the results obtained below by more than a small factor. 
The contribution of the plasma oscillations to the specific heat is then 


. 3 2 
dE vw, _ 45kPomx XE 


dT  3(2mh) ~ (e7—1)° 

where x = E,)/kT. The expression x*e7/(e7—1)? is a monotonic decreasing 
function, having its largest value (for positive x) at x = 0, at which point it is 
unity. It follows then that C, has its largest value when Ey < kT; i.e., 

4mkPo max 


C wes ey 8 


3(2h)° 
Taking the ‘‘number of superconducting electrons’ per unit volume, po/m, 
to be about 5 X10?!, we get: 
Cenex = 5.5 X 10-3 erg/cc. WK. 





C= 





Since measured values (1) of the specific heat in niobium and other super- 
conducting metals are about 104 ergs/cc. °K., we conclude that the con- 
tribution from the plasma oscillations to the specific heat is negligible. 


(2) Electric Field Potential Contribution 
We saw in Section II that the electric field potential & could be expressed 
as 


© = +4, = &)+const. exp(i7K-r+iut)+api 


with K? = ew?/c*?. The contribution of the oscillation p; to the energy of the 
system has been taken into account in III (1). We will now consider the 
contribution of the mode of vibration represented by ,;’. We define the 
energy and momentum of this oscillation in the usual fashion by E = hw and 
P = hK, respectively, and we define the total energy per unit volume of such 
oscillations by 

4e (°E(P)P'dP 4k _ 


to a 7 
(Qah)> Jo eo" —1 15(2xhc)° : 





[19] Ztot = 
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(3) Contribution of the Magnetic Potential 

From equation [16] we see that the magnetic vector potential A can be 
expressed as the sum of the waves grad pi, grad #;’, and the complimentary 
function A’ = Ao exp(iKi-r+iw,t). (¢’, being independent of time, is not a 
wave motion.) The wave represented by A’ can be treated in the same way 
as that represented by #,’. We define the energy and momentum of the A’ 
wave as £; = hw;, and P; = #K, respectively. Equation [17] gives 

Ey? = (c?Py?/e)+£o? 


where Eo = hwo. Hence the total energy of the A’ wave is given by 
4n (°E,(Pi)PidP, 8c” ('° (E,"—Ey)*Ei qe, 


[20] Eitwor = 2. (2ah)® Jo en AT _y = (2ahc)* _ 1 








The factor of two appearing in front of the integral in equation [20] comes 
from the fact that A’ has two possible polarizations, i.e., two independent 
directions for Ay perpendicular to the propagation vector K;. Equation [20] is 
simplified a little by putting uw = E,/kT and uy = E,/kT, and by denoting 


by IJ the integral 
4 
Hes f (1 (A= (ee/s) nid ) u'du ‘ 


Then 


Sr R* ey 
(2xhc)* 


It will be shown below that the specific heat due to the energy of the %,’ 
wave as given by equation [19] is of the same order of magnitude as experi- 
mentally measured values if the dielectric constant is taken to be about 10?°. 
With this e« and with the number of superconducting electrons taken to be 
5X10?!/cc., we get up = 0.05. Hence a very good approximation is made to I 
by taking u» = 0.00, in which case J = 2*/15. We therefore have for the total 
energy of the ®,’ and A’ waves, 





Eiwor = I (uo). 


: 4 5p! 3/274 
[21] Evwt = EvortEitort = “aah . 


Considering the dieiectric constant to be a function of temperature, we get 
for the specific heat contribution from these waves: 


‘ — GE ror _ ( 3/241 =) 
[22] Cc. — dT — aT € +i Ir dT 
where 
> 574 
a= eae = 4.38107" erg/(deg.)‘cc. 


Figure 2 shows a plot of C,*!/aT* versus T calculated from the experimental 
data for niobium (1). This graph shows quite clearly that C,*' does not obey a 
T? law in niobium. 








688 CANADIAN JOURNAL OF PHYSICS. VOL. 33 





Oo i z 3 4 5 t 7 3 TK 


Fic. 2. Variation of C,¢! from a T? law. 


Since we cannot expect to find the functional dependence of the dielectric 
constant on the temperature from this present theory, we must content 
outselves with finding what function will give the closest fit between the 
experimental and theoretical specific heat curves. By a process of trial and 
error, a curve of e versus JT was found which could be used to construct the 
curve in Fig. 2 by the relation 


ce _ 3p P de?” 
ame * +é2 Gr: 





It was found that this empirical curve of e*/? versus T could be described quite 
accurately by the function 


a ee (=) 
{23] e€° = € T xP Tr 


where €9°/? = 16.410" for niobium. 
Hence we can write 


[24] C= a7 Te? (34-22) 

= 1.625X ore 7(3 + a) for niobium. 
Table I shows a comparison of C,*' calculated from [24] with the experimental 
values (1) of C,°'. It should be noted that the experimental values listed for 


T <2.5°K. are extrapolations and not measurements. Hence the large 
differences in the first four entries are not significant. 
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TABLE I 


C;¢!, ergs/°K. cc. 











% difference 














T°. — — 
Experimental Theoretical 

0.5 20 0.00030 

1.0 70 3.23 

1.5 230 97.4 

2.0 700 616 

2.5 1900 2021 +5.9 

3.0 4.56 108 4.74 108 +3.8 

3.5 9.25 X10 9.12 X10? —1.4 

4.0 1.601 X10* 1.528 < 104 —4.8 

4.5 2.458 X 104 2.350 X 104 —4.4 

5.0 3.48 104 3.38 104 —2.9 

5.5 4.69 104 4.64 X10! —1.1 

6.0 6.12 104 6.10 X10* —0.3 

6.5 7.75 X104 7.78 X10¢ +0.4 

7.0 9.59 X10 9.74 X10 +1.5 

7.5 11.7 X10 11.9 X10 +1.7 

8.0 14.2 X10! 14.3 X10‘ +0.7 

8.5 17.0 X10‘ 17.0 X<10¢ 0.0 

8.7 18.1 X10! 18.1 X10 0.0 








IV. CONCLUSIONS 


Demanding that London’s phenomenological equations be consistent with 
the hydrodynamic description of a classical charged fluid has led us to find 
that solutions of the equations of motion contain four independent modes of 
oscillation. One of these, the plasma oscillation, does not contribute appreci- 
ably to the specific heat of the system. The other three oscillations, the one 
electric potential oscillation and the two magnetic vector potential oscillations, 
have been shown to contribute appreciably to the specific heat provided that 
the dielectric constant is, as was suggested by Ginsburg, of the order of 10'°. 

Whether or not this is the correct interpretation of the high specific heat in 
superconducting metals can only be ascertained by an experiment* determining 
the dielectric constant in superconductors. 
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APPENDIX A 
In analogy with the classical Schroedinger field we could start with the 
Lagrangian density: 


[1’] L = uocthe fay wo? -+---(e*—H*) 
8m" p . 8r 


Life: 
we e+v-A)+UG, ®) 
w\C 
*Ginsburg (2) mentions that an experiment has been performed by Galkin and Besugly which 
indicates that the dielectric constant of tin is about 5X10°. The author, however, has been unable 


to find the paper referred to by Ginsburg. 
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where U(p, ©) is some energy density which is to be determined. By indepen- 
dent variation with respect to ¢, p, A, and %, we get the equations of motion: 


(5’] p+V-(pv) = 0, 


. ii? ft ifue)' 1 a aU 
( entaimne en Ga ee aanene 
(7"] viA-S a pv = 0, 
- 49 dU 
(8’] 1 a 
Making the expansions: 
= pPotp1, 
= +4, 
U(p, ®) = U(po, ®o)+A11+Bipit Cipi P1+ $A 2812+ FBopr’, 
where 
ee) 
nee ( db / po, Bo’ Aa \ oa? Jone 
aU a°U 
m(Drae (Boe 
. Op / po, Po Op / po, Po 
7 =) 
es (22 po, Po’ 
and 
ou 
— A,+A2%:+Cipit . 
aU 
ae = Bi+BopitQit+..., 


and assuming that A and ¢ are small such that second order terms in them 
may be neglected, we find that in order to satisfy the zero order equations we 
must have A; = B,; = 0. Hence the first order equations of motion become: 


[5”’] pt py? p——— pov A = 0, 
h? 
[6’’] ea Pi pi = Bopi tC, Fy, 
ad 2 4me 
[7 V A-% pit poV o— ae 5 5 ia 
a = 4 
[8’’] V'e— 4 $; = "(Asti +Cips). 


We can determine A: and C; as follows: take the divergence of [7’’] and demand 
that the Lorentz condition V-A = —(e/c)@ holds. We get 








1- 
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mir 4m 
[7’’a] V *b— 3 @,— = poV *Oo- ans, — =), 
Substituting V°¢ from [5’’] we get 
a ee ee _4ine : 
ae Vo a d; me 2 
Differentiating equation [8’’] with respect to ¢, we get 
< 2: one 4 ° $ 

[8b] Vb: —4 1 = =" (Arb + Cpr). 
Comparison of equations [7’b] and [8b] shows that Az = 0 and C; = —e/m. 
Therefore 

ae = -< pit ... + nigher order terms, 
and 

aU 


— = -* $,+Bopit ... + higher order terms. 
Hence 
U(p, ®) = U(oo, 0) —~ pei t 3 pil < «4 aes cclen theaal 


The absence of any terms in U like (e/m)p:®o or (e/m)po; means that there 
is no Coulomb interaction between the electrons on the average. This means 
essentially that the interaction with the positive lattice charges cancels the 
interaction with the average negative electron charges represented by po and 
). Using these coefficients, the equations of motion become 


[5d] pit povo——— pov A =u, 
> 7 h? 2 é 
[65] +~* ate ©, = Bop. 
| V *A-- SAt ot pve 2? a = 0, 
s al 
[85] Vb — Sot pi = 0. 


We will now demand that London’s first equation is satisfied. Taking the 
gradient of all terms in equation [6b] we have 


h 
~VO+VS, = = v (Bante V ‘n), 
aA ih? 
or ols _p = =e (Bint 4mpo V ‘n) = 0. 


e ) 
Hence (Bits V Aj = 0 
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4m’ p.B 
— pi = C = const. 





or Viet 


The solution of the above equation is: 

4m‘pB: 

= 

If we demand that the spatial average of p; is zero, then it follows immediately 
from [12a] that the constant C must be zero. Hence we have: 

[126] pi = g(t) exp(7Ko-r). 

Proceeding as before we can now construct the solutions to the equations of 
motion. 





a watt exp (iKo* 1) +753 oe 
0 
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THE ABSOLUTE FISSION YIELDS OF TWENTY-EIGHT MASS 
CHAINS IN THE THERMAL NEUTRON FISSION OF U?**! 


By J. A. Petruska,? H. G. THopE, AND R. H. TomLinson 


ABSTRACT 


Twenty-eight absolute fission yields totalling 78% of the heavy and 16% of 
the light fragments have been determined using the mass spectrometer and 
isotope dilution techniques. The precision of the values obtained is in most cases 
better than 2% and the absolute accuracy is estimated to be about 3%. Fine 
structure in the mass-yield curve is discussed in terms of structural preference 
and various chain branching mechanisms. 


INTRODUCTION 


Early radiochemical determinations of fission yields established the asym- 
metric character of the mass~-yield curve for the thermal neutron fission of 
U*%5, However, widely divergent values for the yield of the same mass chain 
were often obtained and, therefore, a smooth curve of best fit was drawn 
through a plot of the mass-yield data such that the sum of the yields would 
total 100%. This work was largely carried out as part of the Manhattan Proj- 
ect and has been summarized by Coryell and Sugarman (4). Thode and Graham 
(23) and Macnamara, Collins, and Thode (16), who reported abnormally high 
yields for Kr*4, Xe'**, and Xe'*4, were the first to establish the existence of 
irregularity or ‘‘fine structure’ in the mass-yield curve. Glendenin (7, 8) 
attempted to explain the high yields of the 133 and 134 mass chains by assum- 
ing that nuclei with 83 neutrons (formed either as primary fragments or as 
secondary products by the emission of a neutron from a primary fragment) 
would have a high probability of ‘‘boiling off’’ a prompt neutron to form the 
more stable 82 neutron configuration. This mechanism was shown to partially 
account for the high yields of masses 133 and 134; but Wiles, Smith, Horsley, 
and Thode (27) have shown that this mechanism fails to completely account 
for the observed fine structure. Wiles (25) and Wiles et al. (27) have therefore 
postulated that, in addition to the Glendenin mechanism, high yields must 
result from an actual structural preference for isotopes with 82 neutrons in the 
primary fission act. Glendenin, Steinberg, Inghram, and Hess (11) reached a 
similar conclusion atter observing a high yield of Mo!°, the complementary 
fragment to the 134 mass chain. 

The more recent radiochemical yields determined by Pappas (20) seem to 
indicate that fine structure in the mass-yield curve is the rule rather than the 
exception. He has also extended the Glendenin hypothesis of neutron “‘boil-off”’ 
to include nuclei with 83, 85, 87, and 89 neutrons in order to explain the yields 
obtained by him for isotopes in the heavy mass region. 

It is difficult to establish from the currently published radiochemical yields 
whether fine structure in the mass-yield curve is widespread. The existence of 
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numerous fission yield values suggests that many are indeed accurate, but it is 
often difficult to distinguish these from the less accurate values. In this regard, 
the 6.62% yield of Xe'** obtained by Katcoff and Rubinson (15) and the 
6.32% yield of Ba!*° obtained by Yaffe, Thode, Merritt, Hawkings, Brown, 
and Bartholomew (29) are perhaps outstanding. The yields of Ba!**, Ba!*°, 
Sr®*, Sr*!, and Mo® have also been determined with considerable precision by 
Reed and Turkevich (22). However, these yields are too few to establish the 
details of the fine structure in the mass-yield curve. 

Relative fission yields obtained with the mass spectrometer are themselves 
usually precise to within a few per cent, but have required normalization to 
some radiochemical yield before their significance in relation to the over-all 
mass-yield curve could be established. In the present work, normalization of 
relative mass spectrometric yields has been achieved without reference to any 
radiochemical yields. Twenty-eight absolute fission yields have been determined 
with the mass spectrometer using isotope dilution techniques. The yields of the 
28 mass chains have been directly related by evaluating the number of atoms 
of a given isotope from the respective isotope dilutions. These values together 
with the measured number of fissions give absolute fission yields. The number 
of fissions was calculated from the change in the B!°/B" ratio in a simulta- 
neously irradiated boron flux monitor and the known cross section for neutron 
fission of U**> relative to the neutron absorption cross section of B!°. Yields 
have been determined for four isotopes of krypton, two of rubidium, and two 
of strontium in the light mass region and four of xenon, three of cesium, two 
of cerium, six of neodymium, and five of samarium in the heavy mass region. 


EXPERIMENTAL 


A piece of natural uranium metal weighing 2.966 gm. (approximate dimen- 
sions 0.8 cm. X 0.5cm. X 0.4 cm.) and a BF; flux monitor were irradiated with 
thermal neutrons in the NRX (Chalk River) reactor. Details of the irradiation, 
separation of the fission products, and evaluation of the neutron flux have been 
described in a previous paper (21). The fission products, after separation from 
the uranium, were deposited on the source filament of a 90° sector type mass 
spectrometer. Chemical separation of the various fission product elements 
from each other was not necessary since each gives positive ions which may be 
distinguished from those of the other elements. Rubidium and cesium yield 
ions at relatively low temperatures of the source filament. At slightly higher 
filament temperatures the mass spectra of these elements disappear and the 
strontium mass spectrum is observed. At higher temperatures lanthanum, 
praseodymium, and neodymium appear as singly charged monoxide ions. As the 
filament temperature is increased still further samarium appears, chiefly as 
singly charged metal ions but also to some extent as metal oxide ions. As the 
filament temperature is further increased, the spectra of these elements grad- 
ually disappear and the spectrum of cerium monoxide ions develops. It has 
not yet been possible to measure the isotopic ratios of the cerium isotopes with- 
out correcting the Ce'*‘ ratio for the presence of Nd!*4. This correction however 
can be obtained since the amount of Nd'‘* may be measured relative to the 
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combined Ce!** and Nd!*4. With the value of the ratio of Nd'** to Nd!*4 
obtained at lower filament temperatures it is therefore possible to subtract the 
contribution of the Nd'*4 and obtain the ratio of Ce!** to the other cerium 
isotopes. The measurement of the ratios of the neodymium and samarium 
isotopes has been described in a previous paper (17). 

In spite of precautions to prevent contamination of the fission products with 
stable isotopes, significant amounts of natural rubidium and strontium were 
found to be present in the determination of the mass spectrometric ratios of 
these fission products. In the case of strontium, the contamination correction 
could be evaluated from the amounts of Sr** and Sr*’? present since these iso- 
topes do not have significant yields in fission. On the other hand, with rubidium 
it was necessary to determine the amount of contamination by means of a 
blank. The procedure for this technique has been described by Wiles, Smith, 
Horsley, and Thode (27). It was further found by means of isotope dilution 
studies that essentially all of the contaminating rubidium and strontium was 
originally present in the uranium rather than in the chemical reagents used 
to separate the fission products. 

Xenon and krypton isotope ratios were not measured as part of this work 
since they had previously been reported by Wanless and Thode (24). In addi- 
tion, isotope dilution was not required for these elements or for cerium since 
the relative yields of each could be normalized through measured yields. With 
regard to the krypton isotopes, the yields of masses 83, 84, and 86 were ob- 
tained relative to the 10.27-year Kr** and the yield of the mass 85 chain was 
determined from Rb**. Similarly, the xenon isotope ratios were determined 
relative to 5.27-day Xe'** and the yield of the 133 mass chain obtained from 
Cs'53,_ The yields of the cerium isotopes were obtained in a similar manner 
from their ratio to 282-day Ce'*‘ and the measured yield of Nd'**, 


RESULTS 


The mass spectrometric data for the elements rubidium, strontium, neo- 
dymium, and samarium are shown in Tables I to IV respectively. The observed 
rubidium and strontium ratios were corrected for contamination by stable 
isotopes originally present in the uranium used. This correction was estimated 
by adding known amounts of uranium from the same original supply to the 
sample and redetermining the above ratios (see Tables I and II). For stron- 
tium the corrected value determined in this way was identical to that obtained 
when one assumed that all the Sr** in the fission product was stable contamina- 
tion and then made the appropriate correction to the Sr** assuming the natural 
abundances of Sr®* and Sr**, The Rb** has been further corrected for the 
fraction of the 85 mass chain which had not yet decayed from Kr**. This cor- 
rection has been made assuming that 0.29 is the branching ratio for the 10.27- 
year Kr*® isomer (1). Since this value requires that 77.5% of the mass 85 chain 
pass directly to rubidium, the ultimate yield of Rb*, after allowing for the 
fraction of Kr** which decays in 2.5 years, is relatively insensitive to the exact 
value of the branching ratio. 
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TABLE I 


THE RATIOS OF THE RUBIDIUM ISOTOPES PRODUCED IN THE THERMAL NEUTRON FISSION OF U% 
TOGETHER WITH ISOTOPE DILUTION DATA 





Mass spectrometric ratio 


Ratio after dilution with natural rubidium 
1. 3.587 X10" atoms per gram of uranium 
2. 4.31310" atoms per gram of uranium 


Atoms of rubidium per gram of uranium* X10" 
if 
2. 
Average 


Ratio after dilution with unirradiated uranium 
(1 part irradiated : 1.81 parts unirradiated) 


Number of atoms of fission rubidium per gram of uranium X10" 


Corrected for natural contamination* 
Corrected for chain branchingt 
Corrected for fast fissiont 


Isotope 
85 87 

0.437+0.001 1.000 
1.207+0.003 1.000 
1.308+0.003 1.000 
7.847 17.96 

7.716 17.66 
7.782+0.07 17.81+0.15 
0.459+0.001 1.000 
7.518 17.70 

9.270 17.70 
9.19+0.08 17.56+0.15 


*Assuming 72.15% Rb® and 27.85% Rb*" in natural rubidium (19). 

tAssuming 0.29 for fraction of mass 85 chain passing through the 10.27-year Kr® isomer (1) 
and 2.55 years between irradiation and separation of the rubidium from the krypton. 

tAssuming 0.8% of each rubidium isotope formed by fast fission of U8, 


TABLE II 


THE RATIOS OF STRONTIUM ISOTOPES PRODUCED IN THE THERMAL NEUTRON FISSION OF U5 
TOGETHER WITH ISOTOPE DILUTION DATA 


86 
Mass spectrometric ratio 
Ratio after dilution with natural strontium* 
1. 2.39510" atoms Sr®8 per gram of uranium — 
2. 2.825 X10" atoms Sr®8 per gram of uranium — 
Atoms of strontium per gram of uranium X10" 
2. — 
Average = 
Ratio after dilution with unirradiated uranium 


(1 part irradiated : 1.81 parts 
unirradiated) 





0.0100+0.0003 0.740+0.007 


0.0241+0.0008 0.858+0.012 


Number of atoms of fission strontium per gram of uranium X10"% 


Corrected for stable contamination* —_— 
Corrected for decay of Sr®°t — 
Corrected for fast fissiont ~- 


Isotope 
88 90 

1.000 

1.376+0.014 1.000 

1.464+0.011 1.000 

2.788 3.768 

2.889 3.905 

2.838+0.050 3.846+0.070 
1.000 

2.516 3.846 

2.516 4.091 

2.496+0.050 4.058+0.070 





*Assuming 82.56% Sr®8 and 9.86% Sr®* in natural strontium (18). 
tAssuming 27.7-year half-life for Sr® (28) and 2.48 years from time of irradiation to time of 


mass spectrometric analysis. 


tAssuming 0.8% of each strontium isotope formed by fast fission of U8, 
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TABLE III 


MASS SPECTROMETRIC AND ISOTOPE DILUTION DATA FOR Nd** pRoDUCED IN THE THERMAL 
NEUTRON FISSION OF U5 


Isotope 
142 143 

Mass spectrometric ratio 0 1.000 
Ratio after dilution with natural neodymium 

1. 1.428 X10'* atoms per gram of uranium 0.664+0.003 1.000 

2. 1.775X10'* atoms per gram of uranium 0.761+0.002 1.000 
Atoms of fission Nd™* per gram of uranium*f 10% 

1 — 4.061 

2. — 4.132 

Average — 4.096 +0 .036 


*Assuming 27.1% Nd'® and 12.2% Nd in natural neodymium (14). 
{Corrected for 0.8% contribution from fast fission of U®*, 


TABLE IV 


MASS SPECTROMETRIC AND ISOTOPE DILUTION DATA FOR Sm"* PRODUCED IN THE THERMAL 
NEUTRON FISSION OF U*% 


Isotope 
148 149 150 
Mass spectrometric ratio 0 1.000 0.255+0.003 


Ratio after dilution with natural samarium 
1. 3.35010" atoms per gram of uranium 0.478+0.004 1.000 
2. 3.963 X10% atoms per gram of uranium 0.516+0.005 1.000 


Atoms of fission Sm"® per gram of uranium*f{ X10" 
1 


5 _- 8.120 — 
2. _ 7.894 —_ 
Average —_ 8.01+0.11 _ 


*Assuming 11.2% Sm* and 13.8% Sm" in natural samarium (13). 

tSince all Sm in fission products is formed by the Sm°(n, -)Sm'® reaction, the total yield of 
Sm? = Sm49+Sm* = 1.955 times the actual number of Sm“® atoms present at time of analysis. 

{Corrected for 0.8% contribution from fast fission of U8, 


The number of atoms of Nd'4* and Sm!*° produced per gram of irradiated 
uranium is evaluated from the mass spectrometric data in Tables III and IV 
respectively. The final values shown for these isotopes as well as those for 
rubidium and strontium isotopes have been reduced by 0.8% to allow for 
contributions from the fast fission of U?**. The 0.8% correction would not be 
expected to apply uniformly to all yields since the mass—yield curve for U** is 
different from that of U?** but, in view of the small magnitude of this correc- 
tion, the error introduced on this account should be small. 

The ratios of the cerium isotopes are shown in Table V. Since the fission 
products were 2.5 years old at the time of analysis, most of the 282-day Ce! 
had decayed to the final product Nd‘. The correction for this loss to the Ce'4 
ratio is therefore large and increases the uncertainty of the yields of the 140 
and 142 mass chains relative to the 144 mass chain. 
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TABLE V 


RELATIVE ABUNDANCES OF THE CERIUM ISOTOPES 
FORMED IN THE THERMAL NEUTRON FISSION OF U?%5 








Mass spectrometric Corrected for 
Isotope ratio decay of Ce!* 
140 1.05 Liz 
142 1.00 1.12 
144 0.101 1.00 





*Assuming 282-day half-life of Ce. 


The relative yields of the isotopes of krypton and xenon taken from the data 
of Wanless and Thode (24) and the neodymium and samarium isotope ratios 
from the data of Melaika, Parker, Petruska, and Tomlinson (17) are shown in 
Table VI. The ratios of the neodymium and samarium isotopes were determined 


TABLE VI 
RELATIVE YIELDS OF XENON, KRYPTON, NEODYMIUM, AND SAMARIUM 
IN THE THERMAL NEUTRON FISSION OF U2 


Krypton (24) Xenon (24) Neodymium (17) Samarium (17) 











Relative Relative Relative Relative’ 
Isotope yield Isotope yield Isotope yield Isotope yield 
83 1.000 131 1.000 143 1.000 147 2.10 
84 1.839 132 1.496 144 0.930 149 1.00 
85 2.387 134 2.750 145 0.665 151 0.398 
86 3.715 136 2.207 146 0.505 152 0.251 
148 0.281 154 0.068 
150 0.110 


for the particular irradiation reported in this work. However, the values given 
in Table VI, which include other determinations, are therefore more reliable. 

Tables I, II, III, and IV show the number of atoms present per gram of 
uranium for each of several isotopes. In order to obtain absolute fission yields 
for these isotopes, it is necessary to know the number of fissions per gram of 
uranium. This determination, which required correction for self-shielding of 
the uranium, has been described in detail in a previous paper (21). The value 
obtained was 7.06X10!* fissions per gram of uranium. The fission yields of 
Rb*®*, Rb§7, Sr88, Sr*°, Nd!48, and Sm" given in Tables VII and VIII are based 
on this value and on the number of fission product atoms per gram of uranium 
for the isotopes given in Tables I, II, III, and IV. The yields of Cs'**, Cs!%5, 
and Cs!*7 were determined for this particular irradiation and the details are 
given in a previous paper (21). The yields of Nd!44, Nd'45, Nd!48, Nd!48, Nd 159, 
Sm!47, Sm?5!, Sm!5?, Sm!54 shown in Table VIII are based on the measured 
absolute yields of Nd'** and Sm!*° together with the relative yields shown in 
Table VI. The yields of Ce!*® and Ce! are obtained from the relative yields 
given in Table V normalized to the yield of the 144 mass chain given by Nd!**. 
The yields of Kr®*, Kr4, and Kr** given in Table VII are based on the relative 
yields given in Table VI and normalized to the measured yield of Rb**. The 
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TABLE VII 


ABSOLUTE FISSION YIELDS OF LIGHT FRAGMENTS 
FORMED IN THE THERMAL NEUTRON FISSION OF U7 





Yield 
Element. ————_—_____ 
This work Literature 
— a 0.15 (4) 
Se — 0.13 (9) 
— = 0.26 (4) 
Kr 0.545 —— 
Kr 1.00 a 
Rb (Kr) 1.30 
Kr 2.02 — 
Rb 2.49 — 
Sr 3.53 
Sr oo 4.78 (22) 
Sr 5.74 = 
Zr oe 5.60 (11)* 
Zr —_— 5.73 (11) 
Ze — 6.12 (11) 
Zr — 6.28 (11) 
Mo 6.35 (11)* 


#4 rbitrarily normalized. 
{Normalized to 6.3% yield for Ba. 


Mass 
chain 








TABLE VIII 


ABSOLUTE FISSION YIELDS OF HEAVY FRAGMENTS 
FORMED IN THE THERMAL NEUTRON FISSION OF U2 
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117-129 
130 
131 
132 
133 
134 
135 
136 
137 
138 
139 
140 
141 


142 
143 
144 


Yield 
Element 
This work Literature 

ae — 1.8(4) 
S oo 2.1 (20) 

Xe 2.92 os 

Xe 4.37 — 

Cs (Xe) 6.59 - 

Xe 8.03 — 

Cs 6.41 — 

Xe 6.44 — 

Cs 6.15 a= 
Ba — 6.4 (22) 

Ce 6.33 —_— 
Ce —- 5.9 (3), 
6.4 (12) 

Ce 6.03 — 

Nd 5.80 a 

Nd 5.39 o 








Yield 
Mass Element. ————---——_—_—— 
chain This work Literature 
96 Zr a 6.26 (11) 
97 Mo — 6.03 (11) 
98 Mo -— 6.00 (11) 
99 Mo c= 6.00 (26)t 
100 Mo _— 6.56 (11) 
101 Mo _ 5.6 (26) 
102 Mo — 4.1 (26) 
103 —_ — 3.1 (4) 
104 — —_— 2.0 (4) 
105 Mo — 1.35 (26) 
106 Ru _ 0.52 (10) 
107-117 — _— 0.4 (4) 
Summation of yields:This work 16.62 
Literature 83.40 
Total 100.0 
Yield 
Mass’ Element — 
chain This work Literature 
145 Nd 3.86 — 
146 Nd 2.93 — 
147 Sm 2.38 _ 
148 Nd 1.63 os 
149 Sm 1.13 —_— 
150 Nd 0.64 — 
151 Sm 0.45 —_ 
152 Sm 0.28; — 
153 — _— 0.15 (4) 
154 Sm 0.077 —_ 
155-162 — _ 0.06 (4) 
Summation of yields: This work 77.84 
Literature 22.8 
Total 100.64 


absolute yields of Xe!*!, Xe'8?, Xe!*4, and Xe!** shown in Table VIII have been 
obtained from the relative yields shown in Table VI normalized to the measured 
yield of Cs!*8. This normalization also requires the ratio of the 133 mass chain 
to the 131 mass chain. The Xe!**/Xe!*! ratio has recently been redetermined 


and is given the value 2.26+0.02 (5). 


In Table VII selected radiochemical yields of mass less than 117 are shown 
in addition to the yields obtained in this work. The mass spectrometric yields 
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of the zirconium and molybdenum isotopes (11) have been arbitrarily nor- 
malized so that the yields in this table total 100%. The yields of the heavy 
fragments shown in Table VIII also include radiochemical values for those 
mass chains where values were not obtained in this work. The data shown in 
Table VII have been plotted in Fig. 1 and the data from Table VIII plotted in 
Fig. 2. Fig. 3 shows Fig. 1 folded over on Fig. 2 in order to facilitate comparison 
of the two. 


DISCUSSION 


The mass spectrometric precision of the relative yields of the isotopes for 
each element is in most cases better than 1%. The yields of Ce!*? and Ce!®, 
however, relative to Ce!*4 may be in error by as much as 5% as a result of the 
large correction for radioactive decay of Ce!44. The yield of Sm!*7 is likewise 
subject to a similar error because of the holdup from its precursor Pm'*7. The 
precision of the values for the number of fission atoms per gram of uranium 
obtained from the duplicate isotope dilutions for each of the elements was 
+1%. It is expected therefore that the relative yield of each fission chain 
obtained for the elements rubidium, strontium, cesium, neodymium, and 
samarium is accurate to within +2%. The ratio of Xe'**/Xe!*! used in the 
normalization of the relative xenon abundances to that of the Cs!** yield has 
been determined with a precision of better than 1%. The xenon yields are 
therefore considered to be well within the +2% limit relating the cesium 
yields to those of the isotopes of the other elements. The krypton abundances 
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Fic. 1. Yields of light fragments formed in the thermal neutron fission of U*. 


© Values obtained in this work. 
X Literature values. 
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Fic. 3. Comparison of yields obtained for light and heavy fragments in the thermal neu- 
tron fission of U%, 


© Values obtained in this work. 
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- - Mass-yield curve for light fragments. 
— Mass-yield curve for heavy fragments. 
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have been normalized through the relative yield of Kr*® and the measured 
yield of Rb*®. Whereas the yield of the mass 85 chain measured from Rb*® is 
not strongly dependent on the branching ratio to the Kr*® isomers, the yields 
of the krypton isotopes relative to that of the mass 85 chain, obtained by cor- 
recting the measured Kr* yield, may be subject to considerable error on this 
account. Independent isotope dilutions of krypton and xenon (2), however, 
confirm the normalization of the krypton isotopes and hence the value of the 
branching ratio. Thus the relative yields of all the mass chains determined in 
this work are expected to be within +2% except the yields of masses 140, 142, 
and 147. In Figs. 1 and 2 all yields determined in this work have been indicated 
by circles and relative precision indicated by a line through the circle where the 
uncertainty exceeds the circle diameter. 

The absolute fission yields are dependent not only on the number of atoms 
of each isotope produced per gram of uranium but also on the value used for the 
total number of fissions. The absolute accuracy of the number of fissions per 
gram of uranium is dependent on several factors which have been discussed 
in a previous paper (21) and which by themselves place a minimum probable 
error of +3.5%. It will be noted, however, that if the value for the number of 
fission atoms per gram of uranium is in error all the measured yields will be 
affected in the same manner. Several radiochemical yields have recently been 
determined with considerable precision, and where these values are available, 
excellent confirmation with the present yields is found. Thus the value of 
6.62% yield for the mass 133 chain obtained by Katcoff and Rubinson (15) 
is in excellent agreement with the value of 6.59% obtained in this work. Also 
Reed and Turkevich (22) have given a value of 6.3% and Yaffe et al. (29) 
6.32% for the yield of Ba!*°, both of which confirm the Ce!*°® yield of 6.33% 
determined in this work. The yield obtained by Yaffe et al., however, should 
probably be increased by some 3-4% to account for a self-shielding correction 
as discussed in a previous paper- (21). The resultant value of about 6.55% is 
probably more reliable for the yield of this mass chain since only 5% accuracy 
is claimed in this work for the yields of the 140 and 142 mass chains. 

A further indication of the accuracy of the absolute yields may be obtained 
by consideration of their sum. In the heavy mass region 77.58% of the fission 
yields have been determined. Only three mass chains with significant yields 
have not been measured. Using the best literature values for these yields a 
total between 100 and 101% is obtained (see Table VIII). It is therefore esti- 
mated that the absolute fission yields reported in Tables VII and VIII are 
accurate to +3%. 

The fission yields determined in the present work indicate only minor fine 
structure other than that in the mass region 132-137. It may be noted in Fig. 1, 
however, that yields of the even masses 84 and 86 appear high in comparison 
with the yields of the odd masses 83, 85, and 87. Also in the heavy mass region 
(Fig. 2) it appears that the yield of mass 144 is relatively high. In view of the 
uncertainty of the 147 yield little significance can be attached to its deviation 
from the smooth curve through other mass yields in this region. Fig. 1, how- 
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ever, shows further fine structure if the fission yields taken from the literature 
are included. 

In order to discuss the values obtained for the various fission yields, it is 
important to review the four processes which may affect the observed cumula- 
tive fission yields. 

1. Structural preference in the fission act for isotopes of a given mass and 

charge. 

2. Neutron emission from the primary fission fragment. 

3. Delayed neutron emission from fission products which have undergone one 

or more 8 decay processes. 

4. Neutron capture by fission products. 

Any structural preference in binary fission must result in equal yields of 
complementary fragments, and hence if both fragments de-excite by the loss 
of the same number of neutrons, the observed fission yields must be identical. 
The asymmetry of the over-all fission yield curve is, of course, a general mani- 
festation of structural preference. The other three processes, which all result in 
chain branching, affect the fission yield resulting from one fission fragment 
but do not affect the yield of the complementary mass chain. Also it may be 
seen that, if these processes lead to an increased yield of a given mass, they must 
always do so at the expense of the yield of some other mass. Neutron absorp- 
tion increases the yield of the mass chain one unit greater than that of the 
neutron absorber. On the other hand, prompt and delayed neutron emission 
increase the yield of the mass chain one unit lower in mass than the neutron 
emitter. The neutron absorption process, however, is dependent on the neutron 
flux during the fission product formation, whereas the neutron emission pro- 
cesses are independent of the flux. These considerations have been reviewed by 
Pappas (20). 

In order to assess the experimental fission yields in terms of the above pro- 
cesses, it is convenient to fold the fission yield curves over each other as shown 
in Fig. 3. Since the sum of the mass numbers of corresponding fission fragments 
together with the number of emitted neutrons must be 236 (U?*5+1 neutron), 
it is necessary to decide which masses should be made to correspond. Arbi- 
trarily, it has been assumed that most fission fragments are de-excited by the 
emission of one prompt neutron, and hence the yields of the light mass chains 
have been plotted in Fig. 3 to correspond to the heavy mass chain such that 
their combined masses total 234. Since the fission fragments do not all lose 
the same number of neutrons in de-excitation, there is some justification for 
this particular matching in spite of the fact that on the average 2.5 neutrons 
are emitted per fission. 

In Fig. 3 it may be seen that the high yield at mass 90 could result from chain 
branching in the mass 91 chain, but any such simple mechanism would merely 
displace the indentation in the mass-yield curve and could not account for the 
apparent displacement of the masses 80-90 in relation to their complementary 
fragments. This increased width of the mass—yield curve of the light fragments 
in relation to that of the heavy fragments could result from a greater tendency 
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towards multiple neutron emission for fragments with masses either less than 
about 95 or greater than about 139. In view of Fraser’s (6) observation that 
30% more neutrons are emitted by fragments in the light mass region, it is 
probable that fragments with mass less than about 95 have greater probability 
of emitting two neutrons. Such neutron emission would not only increase the 
width of the light mass—yield curve in relation to that of the heavy fragments, 
but would also lead to reduced cumulative yields in the mass region where 
transition from one toward two neutron emission occurs. Since the yields of 
masses 91 and 92 appear low in relation to that of mass 90, it is possible that 
many of the primary fission fragments of mass 92 or less de-excite by emission 
of two neutrons, whereas a smaller fraction of the primary fission fragments of 
mass 93 or greater emit more than one neutron. Thus not only would the cumu- 
lative yields of masses 91 and 92 be reduced in relation to those expected on 
the basis of single neutron emission, but so also would the yields of other chains 
of higher mass in this transition region. Correction for this effect could raise 
the yields of several masses by as much as 0.5% fission yield. 

Since some information is available with regard to delayed neutron emitters 
and fission products with large neutron capture cross sections, it is of interest 
to consider what influence these may have on the observed fission yields. If 
fission products are formed at low neutron fluxes and short irradiation times, 
the effects of neutron capture are minimized. Of isotopes in relatively high 
yield, only Xe!*®, Sm14°, and Sm"! have neutron absorption cross sections large 
enough to materially affect the observed fission yields. The corrections re- 
quired by these absorption processes have been discussed in previous papers 
(21, 17) and the yields reported have been corrected for these effects accor- 
dingly. The effects of delayed neutron emission, however, are more difficult to 
assess. Several delayed neutron emitters have been identified by the half-life 
of their emission, but only I!*7 and Br®’ have been given mass assignment with 
any certainty. The relatively high yield of Kr** in relation to Rb*’ and the 
lack of fine structure in the corresponding heavy fragments as shown in Fig. 3 
suggests that the delayed neutron activity of Br’’ could be responsible for the 
observed fine structure at these masses. In a similar manner, other delayed 
neutron emitters must introduce a difference between the yields of the light 
and heavy mass regions. It is suggested that delayed neutron emission by 
As*> may be responsible for the apparently high yield of Kr*4 and low yield of 
Kr*5, Pappas (20) has summarized the extent to which delayed neutron emit- 
ters would affect the mass-yield curve. I'*” and at least two other delayed 
neutron emitters with unidentified mass are sufficiently abundant to alter a 
given mass chain by 0.5% fission yield. Comparison of the light and heavy 
mass regions of the mass-yield curve as shown in Fig. 3 indicates that only 
between masses 132 and 138 could such large alterations occur unless they are 
masked by some form of compensation. From the fission yield data alone, it is 
difficult to predict which of the isotopes in this mass region (in addition to I'*”) 
are the delayed neutron emitters since it appears that other phenomena also 
influence the observed yields. 

It has previously been pointed out (27, 11) that part of the high yields 
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observed at masses 133 and 134 as well as at the complementary masses 100 
and 101 probably result from the preferential formation of nuclides with 82 
neutrons in the primary fission act. Since the yields of masses 133 and 134 are 
considerably greater than those of masses 100 and 101, it is also indicated that 
the yields of masses 133 and 134 are supplemented by chain branching pro- 
cesses. 

Glendenin (7, 8) has postulated that primary fission isotopes reaching the 
83 neutron configuration tend to “boil off” a neutron to reach the more stable 
82 neutron shell structure. Pappas (20) has extended this postulate to include 
nuclides with 85, 87, and 89 neutrons. Such chain branching may supplement 
the yields of masses 133 and 134 but can only do so at the expense of heavier 
mass chains. Since in Fig. 3 the observed yields in this mass region are higher 
than the observed yields of the corresponding light mass fragments, it appears 
that this cannot be the case. It has already been pointed out, however, that 
the observed cumulative yields in the complementary light mass fragments 
may have also been reduced by transition from single to multiple neutron 
emission in this region. The Glendenin—Pappas mechanism, therefore, together 
with delayed neutron emission may well account for part of the enhanced fine 
structure at masses 133 and 134. 

Pappas (20) has attempted to estimate the influence of neutron “‘boil-off”’ 
on the observed cumulative fission yields. He has predicted the presence of 
considerable fine structure in addition to that in the region of mass 134 which 
could extend from mass 136 to mass 144. Although the present data cannot 
definitely exclude this possibility since yields have not been obtained for 
masses 138, 139, and 141, the yields above 140 apparently fall on a reasonably 
smooth curve. Also the accurately determined radiochemical yield at 139 
obtained by Reed and Turkevich (22) and the average of values obtained for 
the yield of mass 141 fit this same smooth curve. In spite of these discrepancies, 
it appears that at least part of the enhanced yields at masses 133 and 134 must 
result from the neutron ‘“‘boil-off’’ mechanism as described by Pappas. The 
quantitative treatment of this, however, must await more nuclear data and 
accurate methods of evaluating the charge distribution of primary fission 
fragments. 
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LETTERS TO THE EDITOR 





Under this heading brief reports of important discoveries in physics may be published. These 
reports should not exceed 600 words and, for any issue, should be submitted not later than six weeks 
previous to the first day of the month of issue. No proof will be sent to the authors. 


a’ Corrections to the Ground State of Helium 


In a recent paper (2) Schrodinger-type wave equations were derived for systems of three or 
more interacting particles. The three-particle case of greatest interest was that of the helium 
atom. Chandrasekhar and Herzberg (1) have found such close agreement between the best 
available experimental value of the energy of the ground state and a new theoretical value 
based on the conventional formulation that it is to be expected that the quantum-theoretical 
corrections of relative order a contribute nothing to the order required. This has been found 
to be the case. Calculations have been carried out for 15 likely terms, which together yielded a 
final correction of the size of 10-2 cm.~!. Since the corrections involved are all proportional to 
Z? they remain small for all atoms of moderate atomic number. Accordingly, it is considered 
unnecessary to report the work in any detail. (The contribution of the Lamb Shift has been 
considered by previous writers.) It is interesting to note that the same result is obtained 
whether one uses equation [20] or equation [31] of Reference (2) below. The methods employed 
in carrying out the calculations were essentially those due to Salpeter (3). 


1, CHANDRASEKHAR, S. and “oy G. Phys. Rev. 98: 1050. 1955. 
2. Ropinson, R. O. A. Can. J. 369. 1955. 


3. SALPETER, E. E. Phys. Rev. 
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